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Abstract 

We consider the normal matrix model with a cubic potential. The 
model is ill-defined, and in order to reguralize it, Elbau and Felder in- 
troduced a model with a cut-off and corresponding system of orthogonal 
polynomials with respect to a varying exponential weight on the cut-off 
region on the complex plane. In the present paper we show how to define 
orthogonal polynomials on a specially chosen system of infinite contours 
on the complex plane, without any cut-off, which satisfy the same recur- 
rence algebraic identity that is asymptotically valid for the orthogonal 
polynomials of Elbau and Felder. 

The main goal of this paper is to develop the Riemann-Hilbert (RH) 
approach to the orthogonal polynomials under consideration and to ob- 
tain their asymptotic behavior on the complex plane as the degree n of 
the polynomial goes to infinity. As the first step in the RH approach, we 
introduce an auxiliary vector equilibrium problem for a pair of measures 
Qui, M2) on the complex plane. We then formulate a 3 x 3 matrix valued 
RH problem for the orthogonal polynomials in hand, and we apply the 
nonlinear steepest descent method of Deift-Zhou to the asymptotic anal- 
ysis of the RH problem. The central steps in our study are a sequence 
of transformations of the RH problem, based on the equilibrium vector 
measure (^1,^12), and the construction of a global parametrix. 

The main result of this paper is a derivation of the large n asymptotics 
of the orthogonal polynomials on the whole complex plane. We prove that 
the distribution of zeros of the orthogonal polynomials converges to the 
measure /xi, the first component of the equilibrium measure. We also 
obtain analytical results for the measure fii relating it to the distribution 
of eigenvalues in the normal matrix model which is uniform in a domain 
bounded by a simple closed curve. 
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1 Introduction 



The normal matrix model is a probability measure on the space ofnxn normal 
matrices M of the form 

-^exp(-nTrV(Af))dM. (1.1) 



Z, 



A typical form for V is 



V(M) = — (MM* - V(M) — V(M*))) , t >0, (1.2) 

to 

where V is a polynomial and V is the polynomial obtained from V by conju- 
gating the coefficients. In this case, the model may alternatively be defined on 
general n x n complex matrices M, see e.g. [33]. In this paper we study in 
particular (ll.l[) - (|1.2l) with a cubic potential 

V(M) = |m 3 , t 3 > 0. (1.3) 

The main feature of the normal matrix model is that the eigenvalues of M 
fill out a bounded two-dimensional domain as n — > oo with a uniform density. 
Wiegmann and Zabrodin [32] showed that if 

oo 

VIM) = y -M k (1.4) 

k=l 

then — Q(to;ti,t2, . . .) is such that 

t Q = -area(fi), t k = --ff ^p, fc= 1,2,3,.... (1.5) 

where denotes the two-dimensional Lebesgue measure in the complex plane. 
The relations (|1.5p characterize the domain f2 by means of its area nto and its 
exterior harmonic moments tk for k > 1 (it is assumed that €= O, and the 
integrals in (|1.5p need to be regularized for /c < 2). An important fact, first 
shown in [32], is that the boundary of Vt as a function of to > (which is 
seen as a time parameter) evolves according to the model of Laplacian growth. 
The Laplacian growth is unstable and singularities such as boundary cusps may 
appear in finite time. See also [27] [30] 133] [34] for related work and surveys. 

However, from a mathematical point of view, the model (fTTT]l - (fTT2|) is not 
well-defined if V is a polynomial of degree > 3, since then the integral 

Z n = J exp(-nTrV(M))dM (1.6) 

diverges, and one cannot normalize the measure (|l.ip to make it a probability 
measure. To make the model well-defined, Elbau and Felder [17l [18] propose 
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to use a cut-off procedure. Instead of considering all normal matrices M, they 
restrict to normal matrices with spectrum in a fixed bounded domain D. The 
integral (|1.6p restricted to all such matrices is convergent and the model is well- 
defined. As it is the case for unitary random matrices, the eigenvalues of M are 
then distributed according to a determinantal point process with a kernel that 
is built out of orthogonal polynomials with respect to the scalar product 

(f,9)D = JJ D f{z)W)e- nV{z) dA(z) (1.7) 

(which depends on n), with 

V(z) = ^(\z\ 2 -V{z)-V{zj), (1.8) 

see [T7]- For each n, we have the sequence (Pk,n)kLo of monic polynomials (i.e., 
Pk,n(z) — z k + ■ ■ ■ ) such that 

(Pk,n j Pj.n) D = ^fc,n<5j,fc, 

and then the correlation kernel for the determinantal point process is 
K n (w, z) = e -f (VW+VW) g MM) , 

k—0 

Elbau and Felder [17l [18] prove that for a polynomial V as in (|1.4[) with ti = 0, 
|*2 1 < 1) and for to small enough it is possible to find a suitable domain D 
such that indeed the eigenvalues of the normal matrix model with cut-off D 
accumulate on a domain Q as n — > oo. The domain Q is characterized by (|1.5|) 
and so in particular evolves according to Laplacian growth in the time parameter 
to- Note that the cut-off approach works fine for to small enough but fails to 
capture important features of the normal matrix model such as the formation 
of cusp singularities at a critical value of to . 

Elbau [T7] also discusses the zeros of the orthogonal polynomials P n ^ n as 
n — > oo. For the cubic potential (| 1 . 3[) and again for to sufficiently small, he 
shows that these zeros accumulate on a starlike set 

£i = [0,a;*]U[0,wx*]U[0,u;V], u = e 2 ™ /3 , (1.9) 

for some explicit value of x* > 0. The set Ei is contained in tt. The limiting 
distribution of zeros is a probability measure on Si satisfying 

J\og\z-(\d^(() = -^-JJjog\z-(\dA(0, zeC\Q. (1.10) 

For general polynomial V and to small enough, Elbau conjectures that the zeros 
of the orthogonal polynomials accumulate on a tree- like set strictly contained in 
Q with a limiting distribution whose logarithmic potentials outside of f2 agrees 
with that of the normalized Lebesgue measure on f2. 
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In this paper we want to analyze the orthogonal polynomials for the normal 
matrix model without making a cut-off. We cannot use the scalar product 



defined on C since the integral (|1.11[) diverges if / and g are polynomials and V 
is a polynomial of degree > 3. Our approach is to replace ()1.11|) by a Hermitian 
form defined on polynomials that is a priori not given by any integral, but that 
should satisfy the relevant algebraic properties of the scalar product ([1.11)1 . 
We define and classify these Hermitian forms. See the next section for precise 
statements. 

It turns out that there is more than one possibility for such a Hermitian form. 
We conjecture that for any polynomial V it is possible to choose the Hermitian 
form in such a way that that the corresponding orthogonal polynomials have the 
same asymptotic behavior as the orthogonal polynomials in the cut-off approach 
of Elbau and Felder. That is, for small values of to the zeros of the polynomials 
accumulate on a tree like set Si with a limiting distribution satisfying (ll.iop . 

We are able to establish this for the cubic case ([1.3)1 and this is the main 
result of the paper. We recover the same set Ei as in Ijl. 9[) and also the domain ft 
that evolves according to the Laplace growth, as we will show. In our approach 
we do not have to restrict to to sufficiently small. We can actually take any to 
up to the critical time io.crtt' F° r this value the endpoints of Si come to the 
boundary of fl. Then three cusps are formed on the boundary and the Laplacian 
growth breaks down. An asymptotic analysis at the critical time would involve 
the Painleve I equation, which is what we see in our model and that we will 
address in a future paper. 

We note that Ameur, Hedenmalm and Makarov [2 [21] do not use a cut- 
off. Instead they consider cases where V(z) tends to +00 as z — > 00 in all 
directions of the complex plane. This approach does not include cases ([1.8)1 
with a polynomial V of degree > 3. 

2 Statement of results 
2.1 Hermitian forms 

We propose to consider orthogonal polynomials with respect to Hermitian forms 
that share the algebraic properties of the scalar product ([1.11)1 with V(z) given 
by ([1.8[1 . The Hermitian form is a sesquilinear form (•, •) defined on the vec- 
tor space of complex polynomials in one variable that satisfies the Hermitian 
condition 



We use the convention that (•, •) is linear in the first argument and conjugate- 
linear in the second. 

To see what kind of condition we want to put on the Hermitian form we look 
at the scalar product (jl.7[) on the cut-off region D. From the complex version 




(1.11) 



(f,g) = (gJ)- 



(2.1) 
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of Green's theorem we obtain for polynomials / and g 

% I f{z)g-{z)e- nV(z) dz = t [[ £ \f(z)gJz)e- nV ^] dA(z) 



2« JdD ' JJd d z 

= t JJ D f(z)'g 7 ^)e- nV ^dA(z)-n JJ^zf(z)gjr)e' nV ^dA( Z ) 

+ n J 'J J(z)V'(z)g(z)e- nV ^dA(z), (2.2) 

which can be written as 

to(f,g') D -n(zf, g ) D + n(f 1 V'g)D= t £ I f {z)W)^ nV {z) dz . (2.3) 

2 « JdD 

The cut-off approach works if the domain D is chosen such that the effect of the 
boundary dD becomes small when n is large. This means that the boundary 
integral in (|2 .3[) is exponentially small compared to the other terms in (|2.3|) in 
the large n limit. 

Inspired by (|2.3[) our idea is to ignore the right-hand side of (|2.3j) and require 
that the Hermitian form satisfies the identity 

to (/, g') - n(zf, g) + n(f, Vg) = (2.4) 

for all polynomials / and g. We call (|2.4I) the structure relation. Combining 
(|2.4| with (|2.1[) we also have the dual structure relation 

to(f, g) ~ n(f, zg) + n(V'f, g) = 0. (2.5) 

We do not impose the condition that the Hermitian form is positive definite, 
and therefore it may not be a scalar product. 

Our first result is a classification of Hermitian forms satisfying (|2.1j) and 
(|2.4| . Sesquilinear forms that satisfy (|2.4j) can be found as follows. We assume 
V is a polynomial and 

d = dcgV -1. 

There are d + 1 directions in the complex plane, given by arg z — 0j , j — 
0, 1, . . . , d, such that 

arg(— V(z)) — > as \z\ — > oo with argz = 9j. 

We choose the directions in such a way that 

2ttj 



7j - U 



1 



j =0,l,...,d+l. (2.6) 



Let Tj for j = 0, . . . , d, be a simple smooth curve in the complex plane, starting 
at infinity in the direction argz = 9j and ending at infinity in the direction 
argz = 9j + i. See Figure [1] for possible contours Tj in case d = 2 and t 3 > 0. 
Then define for polynomials / and g 7 and for j, k = 0, 1, . . . , d, 

*j, k (J,9)= I dz (_ dwf(z)g( W )e-^ wz - nz) - VM \ (2.7) 
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Figure 1: Contours To, Ti, T2 for the case of a cubic potential (d = 2) 

where g is the polynomial obtained from g by conjugating the coefficients. The 
integrals in (|2.7|) are convergent because of our choice of contours. 

It is an easy exercise, based on integration by parts, that the sesquilinear 
forms < I ) j,fe(-, •) indeed satisfy (|2.4p and (|2.5I) . Then so does any linear combina- 
tion. We prove the following. 

Theorem 2.1. For every n € N, to > 0, and polynomial V of degree degV = 
d+ 1, the real vector space of Hermitian forms satisfying (|2.1[) and (|2.4[) is d 2 
dimensional. 

Every such Hermitian form can be uniquely represented as 

d d 

■ f -'-i Y^l r ^Mf.>n (2.8) 

where k is given by (|2.7p and where C = (Cj,fc)^ fe _ is a Hermitian matrix 
having zero row and column sums. 

The proof of Theorem 12.11 is in Section [31 

Now the following obvious question arises. Given a polynomial V, is there a 
Hermitian matrix C = (Cj,k) such that the Hermitian form (|2.8p captures the 
main features of the normal matrix model? We answer this question for the 
cubic case. 



6 



2.2 The cubic potential and rotational symmetry 

We consider from now on the case 

V{z) = |z 3 , i 3 > 0. (2.9) 

Then d = deg V — 1 = 2 and so by Theorem 12.11 the space of Hermitian forms 
(|2.8p is four-dimensional. For the case (|2.9I) it is natural to require an additional 
rotational symmetry 

(f(ujz),g(u;z)) = (f,g), u^e 2 ^ 3 . (2.10) 

The condition (12.10)) corresponds to the fact that V(wz) = V(z) for the cubic 
potential (|2.9I) so that the integral in (jl.lip is invariant under the change of 
variables z H> ojz. Thus (|2.10j) is specific for the cubic potential and will have 
to be modified for other potentials. 

One may verify from the definition (|2.7[) with V given by (|2.9I) that 



$j,fe(/(u;-),fl<w)) = ®j+i,k+i(f,g), 
where the indices are taken modulo 3. Thus 

2 2 2 2 

j=Q k=0 3=0 k=0 

where again the indices are taken modulo 3 so that for example C_i _i = 6*2,2- 
It follows that (|2.8p satisfies the symmetry condition (|2.10j) if and only if 

Cj t k — Cj-i,fe-i, (indices modulo 3). (2-11) 

The condition (|2~TTj) means that C = {Cj t k) is a circulant matrix. As dictated 
by Theorem 12. 1\ we also require that C is Hermitian with zero row and column 
sums. 

The real vector space of circulant Hermitian matrices of size 3x3 with zero 
row and column sums is two-dimensional. A basis is given by the two matrices 



(2.12) 



It turns out that we are able to do asymptotic analysis on the orthogonal poly- 
nomials only if we choose for C a multiple of the second basis matrix in (I2.12[) . 
However, we have no a priori reason to prefer this matrix above the other one, 
or above a linear combination of the two. It is only because of our ability to do 
large n asymptotics that we choose 






(2.13) 



This leads to the following definition of the Hermitian form and the correspond- 
ing orthogonal polynomials. 
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Definition 2.2. Given n, to, £3 > we define the Hcrmitian form (•, •) on the 
vector space of polynomials by 



</>S> = 2^EE £ ^ / d * £ dwf(z)g(w)e-^ wz -^ w3+z ^ (2.14) 




j=Q k=Q 

where 

fe,*)? lfc =o = I 1 -1 I . (2.15) 

Note that the Hermitian form (|2.14|l depends on n even though we do not 
emphasize this in the notation. For each n, we denote by (Pk, n )k=o,i,2,... the 
sequence of monic orthogonal polynomials for the Hermitian form (|2.14p . That 
is, Pk,n{z) — z k + ■ ■ ■ is a polynomial of degree k such that 

(P k , n ,z j )=0, far j = 0,l,...,fc-l. (2.16) 

2.3 Existence of orthogonal polynomials P n n 

The Hermitian form (•.•) may not be positive definite, and therefore the ex- 
istence and uniqueness of the orthogonal polynomials is not guaranteed. This 
will be our next main result. We focus on the diagonal polynomials P n ,n and 
for these polynomials we also determine the limiting behavior of the zeros. 

Theorem 2.3. Let f 3 > and define 

= 8t 



to,crit — ^72" (2-17) 

3 

Then for every f £ (0,to C rit) the orthogonal polynomials P n ^ n for the Hermi- 
tian form (|2.14p exist if n is sufficiently large. In addition, the zeros of P n , n 
accumulate as n — > 00 on the set 

2 

Si = {z e C I z 3 e [0, {x*) 3 }} = \J[0,u) j x*}, (2.18) 

3=0 



where lo = e 2 ^' 1 ^ 3 and 



J- - y/l - 8t *i) ' ■ (2.19) 



Theorem [273] will follow from a strong asymptotic formula for the orthogonal 
polynomials, see Lemma 16. II 



2.4 Limiting zero distribution 

For a polynomial P of degree n with zeros Z\,...,z n in the complex plane, we 
use 

n 

K^) = -E^ ( 2 - 2 °) 

n z — ' 

to denote the normalized zero counting measure. 

The normalized zero counting measures v(P nin ) of the orthogonal polyno- 
mials P n ,n have a limit that we characterize in terms of the solution of a vector 
equilibrium problem from logarithmic potential theory. We use the following 
standard notation. For a measure /z, we define the logarithmic energy 

J(/i) = J J ] og --L-dii(x)d l i(y) (2.21) 

and for two measures \x and we define the mutual energy 

/(/x, v) = fj log —L-dv( X )dv(y). (2.22) 

In addition to the set Si from (|2.18|) we also need 

2 

S 2 = {z e C | z 3 e Mr} = Q [0, -cj j oo). (2.23) 

1=0 

Definition 2.4. Given £0,^3 > 0, we define the energy functional 
£(Mi,M2) = I (pi) +1(^2) -I(m,fJ,2) 



+ U{^r/ im -r') d ^ <2 - 24) 



The vector equilibrium problem is to minimize (|2.24j) among all measures \x\ 
and fi2 such that 

fdfii = l, supp(/ii) C Si, fd(j, 2 =^, supp(^ 2 ) C S 2 , (2.25) 

where S x and S 2 are given by (|2~T8)) and (j2~23|) . 
We prove the following. 

Theorem 2.5. Let t 3 > 0, < t < io.crit arid Zet £* > be given by (|2.19l) . 

T/ien i/iere is a unique minimizer (p*,^) / or ^(^1^2) among all vectors of 
measures (/ii,/i 2 ) satisfying (|2.25j) . J/ to < io.eri* fieri the first component jjb\ 
of the minimizer is the weak limit of the normalized zero counting measures of 
the polynomials P nn as n — > 00. 
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Figure 2: The set Si where the zeros of P n ,7i accumulate and the boundary of 
the domain il for the values to = 1/2 and £3 = 1/4 (left) and for the critical 
values t — 2 and t 3 = 1/4 (right). 

2.5 Two dimensional domain 

We finally make the connection to the domain fi that contains the eigenvalues 
in the normal matrix model, and that is characterized by the relations (|1.5|) . In 
the cubic model that we are considering we are able to construct the domain 
Q in terms of the solution (nl,^) °f the vector equilibrium problem for the 
energy functional (|2.24[) as follows. 

Theorem 2.6. Let t$ > and < to < to,crit- Let /if be the first component of 
the minimizer of the vector equilibrium problem, as described in Theorem \2.5l 
Then the equation 



defines a simple curve d£l that is the boundary of a domain containing Si. 
The domain Q satisfies (11.51) and is such that 



The identity (|2.26[) means that the left-hand side of (|2 . 26[) is the Schwarz 
function of dfl. See Figure [5] for the domain ft 

The proof of Theorem 12.61 is in section @] It follows from an analysis of the 
vector equilibrium problem and an associated three sheeted Riemann surface. 
This Riemann surface will also be important in the proofs of Theorems 12.31 
and 12.51 that are based on a steepest descent analysis of a 3 x 3 matrix valued 
Riemann-Hilbert problem. The connection to the RH problem is explained in 
section [5] and we refer to the discussion in subsection 15.11 for motivation and 




(2.26) 




(2.27) 
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connections with previous work. The actual steepest descent analysis is done 
in section [6] It is quite involved and ultimately leads to a strong asymptotic 
formula for the polynomials P n _ n as n — > oo, see Lemma 16. 1[ from which the 
Theorems 12.31 and 12.51 will follow. 



3 Proof of Theorem [2J] 

We have already seen that any linear combination (|2.8p of the basic forms 
(|2.7I) gives a sesquilinear form that satisfies (12.41) . We first show that any such 
sesquilinear form is characterized by a unique matrix C = (Cj^k) with zero row 
and column sums. 

Lemma 3.1. Any linear combination of the basic forms (|2.7|) can be written 
in the form (|2.8p with a unique matrix C = (Cj,k)j j._ having zero row and 
column sums. 

Proof. The basic forms (|2 . T[) are such that 
</ 



$ J> fc (/' 9) = 0, for every j 



k=0 

and 

d 

®iMf, 9) = 0, for every k. 

3=0 

Thus one may add a constant to a row or to a column of C and obtain the same 
form (|2.8I) . Hence we may assume that C has zero row sums 

a 



Y C j,k = 0, for every j (3.1) 



k=0 

and column sums 



d 

Gj.k — 0, for every k. (3.2) 

3=0 



To show that the conditions (|3.1[) and (|3.2p determine the sesquilinear form 
, we suppose that 

d d 
j=0 fc=0 

where C = (Cj^J has zero row and column sums, and we prove that C = O. 
Recall that V is a polynomial of degree d + 1 which we write as 

V(z) = ^±1^+1 + ^^), degVi < d. (3.4) 



11 



The directions 6j from (|2.6p are such that 

±±tL z d+1 eR-, axgz = 9 i . 
d+1 

We deform the contour Tj in the definition (|2.7[) of &j,k(.f,g) so that it 
consists of the two rays argz = 8j and argz = 0j+\. Then it follows from (|2.7p 
and (|3.3p that for every polynomial /, 

= </, 1) 

= £/ ^/(^^^(C^-C^) / ^e-^^- VW) . (3.5) 
j=0 ^° fc=0 ^ rfc 

where C_x,fc = C^^. Making the substitution z = see ^ in the integral over z in 
(13. 5p and using (|3.4p we obtain 

^ />oo 

V/ /(e^a;)e- ca;d+1 J (a;)da; = O (3.6) 



3=0 



where c = f > 0, and 

to d+1 ' 



^■(x) = e *'eft Vl < e " ,a > £(C 3 -_ llfc - C7 i)fc ) / d«, e -£C e<Wv <<°» (3.7) 



Taking /(x) = a;( d + 1 )*+ r with non-negative integers I and r, we find from 
2M and (Ell 



(3.8) 



d poo 

/ ^ (d+1) ' +r e- ca:£i+ V J W^ = 0, Ud+1 = 
3=0 J " 



Since this holds for every I — 0,1,2..., and since the functions <pj grow at most 
as 0(e C2xd )) as x — > +oo for some C2 > 0, we find from Q3.8P that 

d 

Xl^+i^i^) =°> r = 0, l,...,d 

which in turn implies that 

4> j (x) = 0, j = 0,...,d. (3.9) 

From (|3.7p we know that the function x i— > <fij(x)e T o Vl ^ e ' ' x ^ is a linear combi- 
nation of the entire functions 



x i v 



dwe-^ ei6iwx - v{w) \ k = 0,...,d. 
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The only linear relation between these function is that they add up to zero. 
Because of (I3.7|) and Q3.9P it follows that Cj_i,fc — Cj_k is independent of k. 
Summing over k = 0, 1, . . . , d and using the fact that the matrix C has zero row 
sums, we find that 

Cj-i,k - Cj,k = 0, for j,k = 0,l,...,d, 

i.e., C has constant columns. Since the column sums are zero as well, we get 
C = O as claimed. □ 



Now we give the proof of Theorem 12.11 
Proof of Theorem \2.1l Using the rule that 



f{z)dz = _f(z)dz 
r Jr 



it is easy to obtain from (|2.7[) that 



It follows that (|2.8p with the normalization (|3.ip - (l3.2p satisfies the Hermitian 
condition (|2.1I) if and only if Cj f. = C^j for every j, k, that is, if and only if 

C = C* (C is Hermitian matrix). 

By Lemma 13.11 any sesquilinear form (|2.8I) is characterized by a unique matrix 
C with zero row and column sums. The space of all (d+1) x (d+ 1) Hermitian 
matrices with the zero row/column sum property is isomorphic to the space of 
all d x d Hermitian matrices. Indeed, the restriction of C to the first d rows and 
columns provides an isomorphism. The dimension of this real vector space is 
d 2 . Hence the real dimension of the vector space of Hermitian forms satisfying 
(HP and flUT) is at least d 2 . 

To complete the proof we show that the dimension is at most d 2 . To that 
end, we consider the moments 

l*j,k = (z j ,z k ) 

where (•,•) is a Hermitian form satisfying (|2.4p and (|2.ip . Suppose V(z) — 

d+1 

J2 t z1 with ^ °- Thcn (HI with f( z ) = zi and 9( z ) = zk implies that 
i=i 

d 

ktofij,k-i - nfij+i,k + n X! *i+i/ i j,fe+i = 0. (3.10) 

1=0 

and by (|2~T|) 

Mj.fe = Mfcj- (3-11) 
From p.lOp and p. lip it is easy to see that all moments are determined by 
the moments jij^ with j,k — 0, ...,d— 1. This block of moments yields a 
Hermitian matrix of size d x d, which is determined by d 2 real parameters. 
Therefore the vector space of Hermitian forms satisfying (|2.4I) and (12 . 1[) is at 
most d 2 dimensional. □ 
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4 Proof of Theorem 12.61 



4.1 Vector equilibrium problem 

We begin by analyzing the vector equilibrium problem for the energy functional 
(|2.24p . For every choice of a;* > there is a unique minimizer (/i*, /i^), since Si 
is compact. Both measures are symmetric under 27r/3 rotations. 
Given fi\, the measure [i\ minimizes the functional 

among all measures on £2 with J d^2 — 1/2. This means that is half of the 
balayage of onto £2- Hence /ij always has full support, supp^) = £2, and 
it is characterized by the property that 

2j]og\z-C\di4(Q = Jlog\z-C\di*UQ, zes, (4.1) 

Given /i^ the measure /xj minimizes the functional 

Ml ^ /(Mi) ~ /(/*!,/*$) + yJ (^ |z|3/2 " f **) 

among all probability measures on Si. Thus [i\ is a a minimizer for an energy 
functional with external field, see [10l[29], and it is characterized by the condition 
that there exists a constant I £ R such that 



2 1 log I* -«(<;)- J log|z-CMMl(C) 



-f(^N 3/2 -r 3 ) o iTsT^L^ (4 - 2) 

The support of [i{ consists of a finite union of intervals, in general. 

The number x* is at our disposal. We want to choose it an optimal way as 
described in the following lemma. 

Lemma 4.1. Let < to < to,crit- Then there is a unique value for x* > such 
that 

• fii has full support, i.e., supp(/i^) = Si, and 

• the density of [l\ vanishes at the endpoints uj 3 x* , j — 0, 1,2. 

If to < tQ tCr it the density of n\ vanishes like a square root at the endpoints, while 
for to = to_ cr it it vanishes with an exponent 3/2. 
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To prove Lemma |4~T1 we first assume that we choose x* satisfying the condi- 
tions in the lemma. From that assumption we will find explicit expressions for 
the measures, from which we can indeed check that the conditions are satisfied. 
The proof of Lemma |4~T1 is given in subsection l4.3l The proof will also give that 
x* is given by (|2.19p . 

Theorem 12.61 is proved in subsection 14.41 

4.2 Riemann surface 

The construction of the measures is based on the consideration of a three sheeted 
Riemann surface 1Z with sheets 

TZ 1 = C\E 1 , K 2 = C\(SiUE 2 ), ft 3 = C\£ 2 . (4.3) 

For k = 1,2, the sheet IZk is connected to IZk+i along Efc in the usual crosswise 
manner. Then 72. is a compact Riemann surface of genus zero. 
Define the Cauchy transforms of the measures ^\ and /i 2 

F k (z)=[^Q, z£C\4 fc = l,2. (4.4) 

These functions have the symmetry property 

F k (ujz) =uj 2 F k (z), zeC\£ fc , fc=l,2. 

The conditions (|4.1I) - (|4.2I) together with the symmetry properties, lead to the 
following relations for the Cauchy transforms 

F ll+ (*) + F h .(z) - F 2 (z) = (^=\z\ 1/2 - t 3 \z\ 2 ) , z € S x , 

to VV*3 / (4.5) 

F 2 , + (z) + F 2: -(z)-F 1 (z) = Q, zeS 2 . 

These functions are used to define the ^-functions that play a major role. 
Throughout the paper we use the principal arguments of fractional powers, 
that is, with a branch cut along the negative real axis. 

Definition 4.2. We define 

= t 3 z 2 +t F 1 (z), z€K u (4.6) 



6»w = ^ u/ 2 ^ z;:' : : z , , , ^ (4-7) 



^z 1 / 2 +t (F 2 (z)-F 1 (z)), z eK 2 nS 07 



n 3 



z 1/2 + t (F 2 (z) - z e K 2 n (5i U S 2 ), 



^ 3lj j_i= z i/2_ toF2(z ) ; 2 € ^ 3 n (s, u &). 1 J 

where Sq,Si,S 2 denote the sectors 



S 
Si 

s 2 



— 7r/3 < arg z < 7r/3, 

7r/3 < argz < 7T, (4.9) 

— 7r < argz < — 7r/3. 



15 



Note that by (|4.4p and (|4.6I) we can express the densities of the measure jii* 
in terms of £1. Indeed, by (14.41) we have by the Sokhotskii-Plemelj formula that 



dfi*(x) = — ^ (Fi,-(x) - Fi,+(x))dx, x e supp(^) 
Zm 



which by (|4.6I) leads to 



d»t(x) = ^-r - ( 4 - 10 ) 

ZiriiQ 

(&,+(x) - £i,+ (x))dx, x G supp (//*). (4.11) 



Similarly 

1 / 2Z 1 / 2 

d/i 2(*) = l ±_ 7r + &.+ (*)-&■+(*) ) dz ' ^esupp(^), (4.12) 



with appropriate choice of signs ± and square roots. 

Lemma 4.3. Suppose that supp(/i^) = Ei and that the density of fi* vanishes 
at the endpoints uj j x* for j = 0,1,2. Then 

to < to,crit = rro' (4-13) 

and the following hold. 

(a) The three functions ^1,^2,^3 given by (|4.6[) - (|4.8[) define a meromorphic 
function on the Riemann surface 1Z whose only poles are at the points at 
infinity. 

(b) The functions £j, j — 1, 2, 3 are the three solutions of the cubic equation 

~ hz 2 ^ - (toh + 1 ) 74 + z 3 + A = 0. (4.14) 



where 



A _ l + 20tpt| - 8t 2 4 - (1 - 8^§) 3 / 2 
A ~ 32^1 ' (4 ' 15) 



(c) The algebraic equation (|4.14[) /ias branch points at z = lo j x* , j = 0,1,2, 
where x* is given by (|2.19p and nodes at the values z — u J x, j = 0, 1, 2, 
where 



3 + yj - 8f t 
4ta 



2 

3 > a;*. (4.16) 



TTie values z = uj 3 x* and z = cu^x, j = 0, 1,2, are the only zeros of the 
discriminant of (|4.14p . 
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(d) We have x > x* > with equality only if to = to t0 rit> an d 

£i(<Jx) = b(iJx) = u 2j x, j = 0, 1, 2. (4.17) 

Proof. The conditions (|4.5p imply that the three functions £i, £2, £3 given in 
(|4.6p — (|4.8[1 define a meromorphic function on the Riemann surface. There is a 
double pole at infinity on the first sheet, and a simple pole at the other point at 
infinity. There are no poles at the endpoints uj 3 x* due to the assumption that 
the density of fx\ vanishes at these points. This proves part (a). 

It follows that any symmetric function of the three £j functions is a polyno- 
mial in z. From (14.61) (14.81) and the fact that 



F 1 {z) = z- 1 + 0{z- 4 ), F 2 {z) = \z- 1 +0{z- 5 ' 2 ) (4.18) 
as z — > 00, it is then easy to see that 

'tx{z)+ti{z)+&{z)=Uz 2 , 

+ - - (*o*3 + £) z, (4.19) 

where A is some real constant. Thus ^1,^2,^3 are indeed the solutions of the 



UJ J X 



cubic equation (|4.14[) . The algebraic equation has branch points at z = 
j = 0, 1, 2, and no other finite branch points. This property allows us to deter- 
mine the constant A. 

The discriminant of (14. 14|) with respect to £ is a cubic equation in £ = z 3 
that we calculated with the help of Maple. The result is 

4*!c 3 + (44 + 4 ^3 + m o4 - 8 )C 2 

+ (444 + !8At tl + 12i§£ 3 - 36A + 12i ^ 1 + 4i^ 3 )C - 27 A 2 . (4.20) 

The discriminant (|4.20p has a root for the value £* = {x*) 3 that corresponds 
to the branch points z — u J x* . The fact that there are no other finite branch 
points implies that either is a triple root of (|4.20p or is a simple root and 
(I4.20p has a double root as well. The case of a triple root happens for the values 
(calculated with Maple) 

1 27 3 

to = —o, A= ? . x* = — , (4.21) 

84' 256t|' 4t 3 v ' 

and this is the only possible combination of values with to > (one combination 
with to < also gives rise to a triple root). The values (|4.21l) are the values for 
the critical case. 

We may from now on assume that (|4.20p has one simple root (* > and 
one double root. Since x* is a branch point that connects the first and second 
sheets, we then have that £2(2) — £.i( z ) vanishes as a square root as z — > x* . For 
x > x* , we have that £2 (%) — £1 (x) is real, while according to (14. 10)) we must 
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have ^2,+ (x) — ^i,+ (x) G «K + for x £ [0, x*], since the density of jit* is positive on 
[0, a;*]. This implies that £2 {x) — £1 (x) > for x in some interval (x* ,x*+S) to 
the right of x* . The definitions (|4.6 p -(|4.7 p imply that > &{x) for x > x* 

large enough. Therefore there is a value 

x > x* (4.22) 

such that £i(x) = £2 (20- Then clearly £ = (a;) 3 is the double root of (I4.20[) , and 
z = x , £ = £1 (x) is a node of the spectral curve (I4.14[) . 

Because of the symmetry of the cubic equation (|4.14p in the variables z and 
£, we can interchange the values of z and £, and we get that z — £i(x), £ = x 
is also a node. Thus (£i(x)) 3 is also a double root of (I4.20p . which because of 
the uniqueness of the double root implies (|4.17p for j = 0. Because of 2tt/3 
rotational symmetry we also have (|4.17p for j = 1,2. This proves part (d) of 
the lemma. 

Taking £ = z in (|4.14l) we find from (I4.17P that z — x is a double root of the 
quartic polynomial 

2z 3 ~t 3 z 4 - ( t t 3 + ^)z 2 + A. (4.23) 



h 

Thus x is one of the critical points of the polynomial (|4.23p . which are easily 
calculated to be z = and the two solutions of 2t 3 z 2 — 3z + (t t 3 + j^j = 0. 
Since x > x* > we discard z = 0, and find two possible values for x 



& = ^£EM, Bass ^v^M. (4 . 24) 

This value must be real and so we conclude that 8t$t 2 < 1. This proves f|4. 13[) . 

The corresponding values of A we now find by substituting z — Xj into (|4.23[) 
and equating to 0. The results are 

, 1 + 20M§ - 8^4 - (1 - 8t t§) 3/2 



32i 3 



An 



20V| - 8t 2 tj + (1 - 8t t 2 3 ) 3 / 2 
32t| 



For the value A — Aj, the simple root of (|4.20p comes out as (* — (x*) 3 with 

x{ = J- (l - yjl-8t t£) 7 , x\ = J- (l + ^1-8*0*1) 7 . (4-25) 

which is again calculated with Maple. From (|4.24p and (|4.25p it follows that 

3 ~ 3 ^ 

X,< — <X!, X a > — >X 2 . 

The inequality (|4.17l) is only satisfied in the first case, and so x* — x*, A = 
A%, and x = x\. This proves the formulas (I2.19p . (|4.15l) . and (I4.16p . which 
establishes the parts (b) and (c) and completes the proof of the lemma. □ 
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Figure 3: Density of n\ on [0, x*] for the values to = 1/4 and ts — 1/2 (left) and 
for the critical values to = 1/2 and £3 = 1/2 (right). In the non-critical case the 
density behaves as c(x* — x) 1 / 2 as x — > x* , while in the critical case it vanishes 
like c(x* — x) 3 / 2 . 

4.3 Proof of Lemma 14.11 

We now prove Lemma |4. II by reversing the arguments given above. 

For to < to,crit, we start from the spectral curve (|4.14[) with the value of A 
as in (|4.15p . The spectral curve defines a Riemann surface with sheet structure 
as in (14.3[) . This defines in particular the value of x* . There is one solution of 
(|4.14|) that satisfies 

(,i(z) = Hz 2 + toz^ 1 + 0{z~ A ) asz^oo 

and £1 is defined on the first sheet. The analytic continuations onto the second 
and third sheets are then denoted by £2 and £3, respectively. 

Then we define a measure /i* on £j for j — 1,2, by the formulas (14.101) and 
f|4. 12|) . All parts of £1 and £2 are oriented away from the origin. The vector 
of measures (ji* , is the minimizer of the vector equilibrium problem, and fi* 
satisfies the conditions of Lemma 14. 11 since it is positive on [0, x*) and vanishes 
at x* , see also Figure [3] for the density of n\ for a non-critical value to and 
for the critical value to, rit- in a non-critical case the density of vanishes 
as a square root at x*, while in the critical case it vanishes like (x* — x) 3 / 2 as 
x — > x* — . See Figure |3] for plots of the density of fx* in the non-critical and 
critical cases. 
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4.4 Proof of Theorem [2^1 

The spectral curve (|4.14[) with A given by (|4.15p has genus zero. The curve has 
a remarkable rational parametrization 

z = h(w) = rw + aw~ 2 , £ = = aw 2 + ra -1 (4.26) 

with 

1-VT^SI a _lz£EM. (4.27, 



2t 3 4t 3 

This can easily be checked by plugging (|4.26|) - (|4.27j) into (|4.14|) . The calcula- 
tions show that it only works if A has the value as in (|4.15|) . 

Assuming to < fo )C Ht we see that T and a are both real and positive. We 
view z = h(w) as a mapping from the w-plane to the z-plane. The mapping is 
conformal around infinity. For large enough p we have that the circle \w\ = p 
is mapped to a simple closed curve. This will continue to be the case if we 
decrease p until the circle |to| = p contains a critical point, that is, a solution of 

h'(w) = r- 2aw~ 3 = 0. 

The critical points are on the circle with radius 

'2a\ 1/3 / I \ 1/6 

' = ( 1 - < 1. 



Peril = [ — J = [ J - \/ 1 - 'Vn'H 

Since p cr a < 1, we see that h(\w\ = 1) is therefore a simple closed curve 
which is smooth for to < to.cHt and has three cusp points for the critical value 
to,crit- Let ft be the bounded domain that is enclosed by h(\w\ = 1). 

For z € <9fi, we have z — h(w) for some w with \w\ = 1. Then £i(z) = 
h(l/w), see (|4.26p . and so since 1/w = w, and since the coefficient of h are real 
numbers 

£i (2) = h(w) = = z, z G 90. (4.28) 

This shows that (|2.26l) indeed defines a curve dft, by virtue of (|4.6I) and (|4.4I) . 
We compute the area of by means of Green's formula 

area(O) = // dA — — <L zdz = — <£ £x(z)dz. 
J Jo. 2 *Jdn 2% ha 

Here we used (|4.28[) . Now recall that ^1 is analytic in the exterior of fl with 

£i(z) =< 3 z 2 + < z _1 + O(z _4 ), as z -> 00, (4.29) 

see (|4.6p and (|4. 18[) . Thus we can move the contour to infinity and by doing so 
we only pick up a residue contribution at 00. This proves the formula for i in 

(E3). 
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The exterior harmonic moments (| 1 . 5[) are computed in a similar way. We 
have by Green's formula applied to the exterior domain 

--[[ ^ = ^-<f 4*=^/ 

n JJc\n z 2m Jan z 2m Jan z 

By contour deformation and using (I4.29|) we find that this is t 3 for k = 3 and 
for other k > 2. Thus ([13]) follows. 

We finally prove (|2.27[) with yet again similar arguments. For z e C \ O we 
have by Green's formula and (|4.28[) 

nJJn z ~C ^Jan z ~C 2m Jan z ~C 

Moving the contour to infinity we pick up a residue contribution at £ = z, which 
is £1(2), and at £ = 00, which is — ^z 2 because of (|4.29[) . In total we find that 

which is equal to t Q Fi(z) by P~o) . Then dividing by t we obtain (j2T2"71) . This 
completes the proof of Theorem 12. 61 

5 Preliminary steps towards the proofs of The- 
orems 12.31 and 12.51 

5.1 Discussion 

The orthogonality induced by the Hermitian form (|2.14j) is similar to the biorthog- 
onality that plays a role in the two- matrix model, see [4] . The two-matrix model 
is a model for two random matrices with two potentials V and W on R with 
sufficient increase at ±00. The biorthogonal polynomials in this model are two 
sequences (pj, n (x))j, degp JiTl = j and (qk,n)k, degg fe ,n = k of monic polynomials 
satisfying 

/OO /"OO 
/ p 3 A*)<lkAy)e~ n{nx)+W{y) ~ TXv) dxdy = h k 6 3 , k (5.1) 

where r is the coupling constant. Comparing (I5.1[) with (|2.14p one sees that 
f|2 . 14[) is like the biorthogonality in the two matrix model with equal cubic 
potentials V(x) = W(x) — —-^f-x 3 and coupling constant r = — j-- The main 
difference is that integrals in (|5.ip are over the real line, while integrals in (|2.14l) 
are over combinations of the contours 1^, j — 0, 1, 2. 

This docs not play a role on a formal level. The main algebraic properties 
that are known for the biorthogonal polynomials, see e.g. [H [5j |6], also hold 
for the orthogonal polynomials with respect to the Hermitian form (|2.14[) . This 
includes the existence of differential and difference equations, and an integrable 
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structure of r-functions and Toda equations. However, analytic and asymptotic 
properties depend crucially on the precise definition of contours. For example, 
it is known that the biorthogonal polynomials characterized by (|5.1[) have real 
and simple zeros [16] . which is not the case for the orthogonal polynomials for 

Recently [HI [T^l HE] , the biorthogonal polynomials p n ^ n from (|5.1[) with an 
even polynomial V and a quartic potential W(y) = \y A + ^y 2 were success- 
fully analyzed in the large n limit by means of a steepest descent analysis of a 
Riemann-Hilbert problem. The Riemann-Hilbert problem was obtained earlier 
in [23] from a reformulation of the biorthogonality (I5.1[) as multiple orthogo- 
nality, since for multiple orthogonal polynomials a Riemann-Hilbert problem is 
known [31 as a generalization of the well-known Riemann-Hilbert problem for 
orthogonal polynomials [20]. Because of the formal similarity the same setup 
works in the case of orthogonal polynomials for the Hermitian form (|2.14[) . We 
can reformulate the orthogonality as multiple orthogonality and it leads to a 
3x3 matrix valued Riemann-Hilbert problem with jumps on the contours Tj. 
In order to prepare for the asymptotic analysis we first adjust the contours in a 
suitable way. This will be done in this section. 

The systematic asymptotic analysis of Riemann-Hilbert problems is due to 
Deift and Zhou [I3J who developed their steepest descent analysis first in the 
context of large time asymptotics of integrable systems. It was applied to or- 
thogonal polynomials and random matrices in [5J [TTJ [T3] . These papers also 
emphasized the use of equilibrium measures in the asymptotic analysis. The 
Riemann-Hilbert problem for orthogonal polynomials [20 a is of size 2x2. Ex- 
tensions of the steepest descent analysis to larger size Riemann-Hilbert problems 
were first discussed in [9] [26] . In the next section we will build on these and 
later works (see [33] for an overview) and apply the Deift/Zhou steepest de- 
scent analysis to the 3x3 matrix RH problem 15.41 A crucial role is played 
by the minimizer (a**, ^2) °f the vector equilibrium problem. This is inspired 
by [7] Q31 [TS] where the steepest descent analysis also depended crucially on a 
vector equilibrium problem. 

5.2 Multiple orthogonality and Airy functions 

The orthogonality with respect to the Hermitian form (|2.14[) is very similar to 
the biorthogonality that plays a role in the two-matrix model. We will use ideas 
that were developed for the asymptotic analysis of the two-matrix model with 
a quartic potential |14[ 115) and apply these to the orthogonal polynomials for 
the Hermitian form (|2.14j) . 

First of all we identify the orthogonal polynomials as multiple orthogonal 
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polynomials. We define for j = 0, 1, 2 and n € N the entire functions 

S&'-^+^dw, (5.2) 



„(z) = -L ej,k [_ we-% {wz -^ w3+z3 »dw. (5.3) 

Lemma 5.1. The monic orthogonal polynomial P n , n is characterized by the 
properties 

2 

V / P n ,n(z)z k W Q .j,n(z)dz = 0, fc = 0, 1, . . . , \% ] - 1, 

7 = 

(5-4) 

V / P„, n (z)z*u;i, J >(*)<fe = 0, A: = 0, 1, ... , Lf J — 1- 
J=o Jr i 

Proof. This follows as in |24j . but for convenience to the reader we present the 
argument. Recall that P n , n satisfies the orthogonality conditions (j2.16[) with 
k = n. 

We prove that for every k, 

(z k P n , n (z),z j ) = 0, j = 0,...,n-2A;-l. (5.5) 

For k — 0, the condition (15 ,5|) reduces to the orthogonality condition (|2.16p . 

Assume that (|5.5[) holds for certain fc > 0. The structure relation (|2.4I) with 
cubic potential (|2.9[) gives 

n(z/,g>=< (/ ! .9')+^3(/^ 2 5)- (5-6) 

Taking /(z) = z k P n . n (z) and = z^ we find that both terms in the right- 
hand side of (|5.6[) vanish for j + 2 < n — 2k — 1, because of the induction 
hypothesis. Then the left-hand side vanishes as well, and this gives (|5.5|) with 
k + 1. Thus (|5.5j) follows by induction. 
Taking J = and j = 1 in (I5.5[) , we find 

(z k p n , n , i) = o, fc=o,...,rn-i, 

(z fc P„^,z)=0, fc = 0,...,LfJ -1, 

and these conditions are equal to the conditions (|5.4[) because of the represen- 
tation (|2 . 14[) of the Hcrmitian form and the definition (|5.2l) - (|5.3p of the weight 
functions. 

The conditions ()5.4p in fact characterize the monic orthogonal polynomial 
since we can similarly show that the conditions ()5.7j) are in fact equivalent to the 
orthogonality conditions (|2.16l) (again by using the structure relation (|5.6[1 ). □ 
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The functions (I5.2[) and (I5.3|) can be expressed in terms of Airy functions 
and their derivatives. The classical Airy differential equation y" = zy has the 
three solutions 



1 f . 1 



y (z) = M(z) = — e -* s - zs ds = -— e 3S ~ zs ds 



2tti ./r 2ni J r„ 



i,3. 



1 „3 



yi(«)=wAi(wz) = — / e^'-^ds = -— / e t s ~ (5.8) 

27TI 7 Fl 27TI Jp 2 

1/2(2) = Ai(w 2 z) = — I ei s3 ~ zs ds = : / e^ e3 ~ zs ds, 

that are related by the identity yo + y\ + 2/2=0. Then we get from (15.21) and 

1/3 



after the change of variables w = (^^^j 



s, 



Wo,o,n(z) = d n (y 2 (c n z) - yi(c n z))e 3t o ' 



wo,i,n[z) = d n (yi(c n z) - y Q (c n z))e 3t o z , ( 5 -9) 



£3,3 



with constants 



Wo,2,n(z) = d n (y (c n z) - y2(c n z))e 3t o 



r, 2 / 3 / i \ 1/3 



Similarly from (|5.3p and (|5.8|l 



*3 ,3 



wi,o,n(«) = ~d n (y 2 (c„z) - yi(c„z))e 3 *o 2 , 



= "^(yU^z) - ^(c„z))e^o z3 , (5.11) 



wi,2,n(z) = -d n (yg(c„z) - y' 2 (c n z))e 3t o 



la. .3 



5.3 Riemann-Hilbert problem 

Observe that there is no complex conjugation in the multiple orthogonality 
conditions (|5.4[) . Following [3T] we then find a characterization of P n , n in terms 
of a 3 x 3 matrix valued Riemann-Hilbert (RH) problem. We assume that the 
contours Tj are disjoint, and have the orientation as shown in Figure [TJ The 
orientation induces a + and — side on Tj, where the + side (— side) is on the 
left (right) while traversing the contour according to its orientation. 

RH problem 5.2. Let T = [f j=Q T r We look for Y : C \ T C 3x3 satisfying 

• Y is analytic in C \ T, 

• Y + {z) = Y-{z)Jy{z) for z G V with jump matrix Jy given by 

/ 1 W ,j, n Wl,j,n\ 
J Y = 1 on Tj, (5.12) 
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(here Y + and Y- denote the limiting values of Y, when approaching V on 
the + side and — side, respectively), 

• Y(z) = (7 + 0(l/z))diag(z",z-r»/2l ;Z -L«/2j) a s z ^ oo. 

The jump condition f!5 . 12[) has to be adjusted in the case of overlapping 
contours. 

Provided that the orthogonal polynomial P n ,n uniquely exists, the RH prob- 
lem 15.21 has a unique solution, see [21] ■ The first row of Y is given by 



Y 1A {z) =P n , n (z) 

/7T1 — ^ 



j=0 
2 



/.in — ' 



2ui 



r, 



Pn,n (^)^0,j,n 

(x) 



P„ lV ,(x)wi t j tn (x) 



dx, 



dx. 



The remaining two rows of Y are built in a similar way out of certain polynomials 
of degrees n — 1 (one of which is proportional to P„_i in , if this orthogonal 
polynomial exists.) 

5.4 Deformation of contours 

In order to prepare for the steepest descent analysis of the RH problem, we 
first adjust the contours Tj. Because the weig hts ([53]) . (|5~TT1) in the multiple 
orthogonality conditions (|5.4|) are entire functions, we have the freedom to make 
arbitrary deformation for each contour Tj , as long as the contours start and end 
at the same asymptotic angles at infinity. We are going to deform the contours 
such that their union contains the set Si = U^_ o [0, ^x*] where the zeros of 
the polynomials are going to accumulate. The deformed contours are shown in 
Figure |U 

We recall that x > x* , see (|4.16|) . We deform To to contain the intervals 
[0,w 2 ir] and [0,wa:] and we continue Tq with an unbounded contour from lux to 
infinity at asymptotic angle 57r/12, and its mirror image in the real axis which 
goes from ui 2 x to infinity at angle — 57r/12. The number 57r/12 could be replaced 
with any number between 7r/3 and 7r/2. We obtain Ti and T2 by rotating Tq 
over angles 2n/3 and — 27r/3, respectively. 

It will be convenient to give new names to the unbounded contours. We 
define 

C+ = (Ti n So) \ [0, x], Co = (T 2 n So) \ [0, x], 

c+ = (r 2 nS 1 )\[o,u)x], c 1 - = (r n5i)\[o,ws|, (5.13) 
c+ = (To n s 2 ) \ [o,u; 2 x}, c 2 - = (r x n s 2 ) \ [o, w 2 x], 



and 



Ci = CfuC7, J =0,1, 2. (5.14) 
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Figure 4: Deformed contours Tq, Ti, T2, All contours are oriented away from 0. 

Then, see also Figure HI 

T = C 2 + U [0,u 2 x] U [0,ux] U Cf , 

ri = C + U [0,x] U [0,uj 2 x] UC^T, (5.15) 
T 2 = C*^ U [0, ujx] U [0, i]UC ", 

and 

2 

r= |J ([oV^UCj). (5.16) 

i=0 

We redefine the orientation on the contours such that all parts are from now on 
oriented away from and towards 00 , as shown in Figure @] We have a freedom 
in the choice of the precise location of the contours Cj and we make use of this 
freedom later on. 

After this deformation of contours and partial reversion of orientation we 
find new expressions for the weights. Indeed, for x £ [0, x], we find by (|5.9|) 

- l"0,l,n(l) + W0,2,n(x) 

= d n (-(j/i(c„x) - y (c n x)) + (y (c n x) - y 2 (c„a;))) eW 1 
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Using the identity yo + y\ + 1/2 — and j/q = Ai, this reduces to 

nt 3 3 

- Wo,i,n(%) + w 0<2 , n (x) = 3d n Ai(c n x)e 3t o x . (5.17) 
Similarly we get from (|5.11l) 

- Wl,l,n{x) + Wl,2,n{x) = -3d 2 ki (c n x)e ~^ x . (5.18) 

It is thanks to our choice (|2.13j) for the matrix C, that we obtain the Airy 
function Ai and its derivative in (|5.17|) - (|5.18|) . and not some other solution of 
the Airy differential equation. The Airy function Ai is a special solution of 
y" = zy because of its asymptotic behavior, see [TJ formulas 10.4.59, 10.4.61], 

Ai(z) = ^e-§* 3/2 (l + 0(z- 3 / 2 )), 

Ai'(zH-f^e-t- V2 (l + 0(z- 3 / 2 )), 

as z — > 00, — 7r < argz < 7r, which is decaying for z going to infinity on 
the positive real axis. Any other linearly independent solution of the Airy 
differential equation increases on the positive real axis. We use (|5.19[) in the 
steepest descent analysis of the RH problem, and this explains why we chose 
the matrix C as we did in (|2.13|) . 

We now redefine the weights on the new contour T by combining the con- 
tributions of overlapping contours. We also rescale the weights by dropping the 
irrelevant constant prefactor 3d n in f|5 . 1 7[) and —3c? 2 in (|5.18j) . 

Definition 5.3. Let T = Uj=o (IP> U Of) ^ e as m FigvrrelU We define two 
functions wo tTi and wi <n on T by first defining them in the sector So ■ | argz| < 
7r/3 by putting 




Ai(c„a;)e 3t o x , 



for 2; G [0,£], (5.20) 
for zeC£, (5.21) 



= 3(2/0 (c n z) -yi(c n z))e 3t o z , 

71*3 3 

7xt 3 3 

= \{ya{c n z) - y2(c n z))e^ z , 

n,t 3 3 

and then by extending to the parts of T in the other sectors by the property 



for zeC(7, (5.22) 



Wi, n (UJz) = 0JWi >n (z), 
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With the definition (15T2"0|) we have by that for x e (0,x], 

-<•> = (-5 (aV - M) f 1 + °<""» • 

(5.24) 

as n — > oo. The exponential part in the asymptotic behavior (j5.24[) of the weight 
functions is reflected in the term 



1 



^ | |3/2_*3 3 \ d {) 



to J 

that appears in the energy functional (|2.24l) . 

5.5 Riemann-Hilbert problem after deformation 

After the contour deformation wc find that the multiple orthogonality conditions 
(|5.4p satisfied by the monic orthogonal polynomial can also be stated as 



(5.25) 



r 



P n ^{z)z k w ,n{z)dz = 0, k = 0, . . . , m - 1, 
P n , n (z)z k w hn (z)dz = 0, k = 0, . . . , Lf J - 1. 
The corresponding RH problem is as follows. 



RH problem 5.4. We look for Y : C \ T -> C JX J satisfying 

• y : C \ T -> C 3x3 is analytic, 

(1 W 0t n{z) Wi tn (z)\ 
1 for z e r 

1 J 

• Y(z) = (7 + 0(l/z))diag(z",z-r™/2l jZ -L«/2j^ as z ^ 



The RH problems 15.21 and 15.41 are equivalent. Also RH problem 15.41 has a 
solution if and only if the monic multiple orthogonal polynomial uniquely exists 
and in that case one has 

Y n (z) =P„, n (z), (5.26) 

as before. 

6 Riemann-Hilbert steepest descent analysis and 
the proofs of Theorems 12.31 and 12.51 



In this section we perform an asymptotic analysis of the RH problem 15.41 For 
convenience we assume that n is even, although this is not essential. The asymp- 
totic analysis will lead to the proof of the following result. 
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Lemma 6.1. Let t 3 > and < to < /o.crit- Then for n large enough the 
orthogonal polynomial P n ^ n exists and satisfies 

P njl (z) = Af M (z)e n9l ^(l + 0(l/n)), as n4 oo (6.1) 

uniformly for z in compact subsets of C \ Si. .Here <?i is defined by 



gi(z) = J log(z-s)df4(s), zeC\Si, 

see oZso (|6.16|) below, and M\.\ is an analytic function with no zeros in C\Si. 

We write M^i since this function will arise as the (1,1) entry of a matrix- 
valued function M, see Section! 



6.1 First transformation 

We start from the RH problem 15.41 for Y with n even. The first transformation 
Y H> X will have a different form in each of the three sectors (|4.9[) . 

We recall the Airy functions yj, j — 0, 1, 2 from (15. 8[) . We also introduce 

yi,{z) = 27ri(w 2 Ai(wz) — a; Ai(w 2 z)), 

2/4(2) = w«/ 3 (wz), (6.2) 
y 5 (z) = uj 2 y 3 (aj 2 z). 

These solutions of the Airy differential equation are chosen such that the Wron- 
skians W{y ,y k ) = y 3 y' k - y' 3 y k satisfy 

W( yj ,y j+3 ) = 1, j — 0, 1, 2, 
W(y 3 ,y 4 ) = W(y 3 ,y 5 ) = W( Vi , y 5 ) = 0. 



(6.3) 



We obviously have from (|6.2I) that 2/3(0) — 0. Hence, by the uniqueness of the 
solution of the initial value problem for the Airy equation at z = 0, we obtain 
that 2/3 (wz) = 2/3(2). This implies that 3/4(2) = ujy 3 {z) and 2/5(2) = ui 2 y 3 {z). 
We also recall the Wronskian relations [TJ formulas 10.4.11-10.4.13] 

W(y ,yi) = W(y 1 ,y 2 ) = W(y 2 ,y ) 1 



2ni 



(6.4) 



Definition 6.2. We define 



A(z) = Y(z) x < 



A 











2/3 ( C n Z ) 


-2/o( C n Z ) 




-2/3 (c n ^) 


2/0 (c„ 2) 










6 


2/5 ( c » z ) 


-2/ 2 ( C n^) 




-2/5 (c„2) 


2/2 (C„2) 


A 











2/4( c n2) 


-2/i (c n 2) 




-J/4 (C„2) 


2/1 (c«2) 



2 e So, 



2 e Si, 



zeS 2 



(6.5) 
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s 2 




Figure 5: Contour Tx for the RH problem for X and the sectors Sq, S\ and S 2 - 



where c„ is the constant given by (j5 . 10[) 

Since the formula (|6.5I) is different in the three sectors, we will have that X is 
discontinuous on the boundary of the sectors, which are the rays arg z = ±7r/3, it 
that form the contour £2 given in (|2.23j) . Thus X is defined and analytic in 
C\ Tx where Fx = ru £2, see Figured! with a jump X + = X-Jx on Tx- The 
jump matrix on the intervals (0,u} J x] simplifies to 



as a consequence of the definitions (16.5[) - (|6.6|) and the Wronskian relation f|6 . 3[) . 



and 




(6.6) 




x 




3 
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The jump matrices on the new contours £2 are piecewise constant, since 
they contain Wronskians W{yi,%)j) of solutions of the Airy equation. Indeed, 
we have by (|6.5p and 



'10 \ I 

J x( z )=(° Vo(c n z) Voicnz)] y' b {c n z) -y' 2 (c n z) 
^0 y 3 {c n z) y' 3 (c n z)J \0 -y 5 (c n z) y2(c n z) 

'10 

W(y ,y 5 ) W(y 2 ,yo) 
v W(y 3 ,y 5 ) W(y 2 ,y 3 ) / 

'10 

w2 ~2^7 I ' for argz = tt/3. 



,0 u 



and then by (|6.6 




Jx(z) = 1 , for arg z = tt/3. 



The jump matrix Jx turns out to be exactly the same for argz = — 7r/3 and 
arg 2: = 7T. 

On the remaining parts of I\x; (the unbounded contours Cj ) , the jump matrix 
takes the form 

1 ae 3t o /3e 3t " ' 
Jx(z) = ( 1 
1 

for certain constants a and j3, that again come from Wronskians of solutions of 
the Airy equation. Let 

2/6 (z) = ^2(2) - 2/1(2)), 2/7(2) = w!/ 6 (wz), y 8 (z) = w 2 y 6 (a; 2 z). 
Then the constants turn out to be 



'a = W(-i/8,ite) = | + |n/3, 
^ = 27rW(-2/ 8 ,2/o) = -|, 

a = W(y 7 ,y 3 ) = \- i«V3, 



on C^, 



2 6 
1 
3 > 



t onC„ , 

/3 = 27riW(j/ 7) y )- 1 



with the same expressions on the other unbounded parts. 
Thus we obtain the following RH problem for X. 

RH problem 6.3. The matrix-valued function X defined by (|6.5[) - (|6.6[) is the 
solution of the following RH problem. 

• X is analytic in C \ Tx , 
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• X + = X- Jx on Tx with 



Jx(z) = < 




As 







oo we have 



1 





1 



X(z) = {I + 0{z^))A{z) 
' (z 11 



z ~ n / 2 e 3ta ^ / 



.3/2 



,3/2 



x < 



0™ 








z~"/ 2 e 3t o\A3' 

V o o 



,3/2 



•»/ 2 e 3*0\/*3 J 

o \ 


, z 3/2 

-n/2g3t 0v /tj / 



with 



'10 
z 1 '* 
,0 



x < 



Q 





i 








i 




V2 


75 












^2 




1 


z 


s/2 


%/2 








i 


i 






1 


i 




V2 



for z £ \JA0,Ujx], 



for z € E 2 , 



for z € U, Cf 



(6.7) 



2 £ Sn 



(6.8) 



for z £ Sn 



~J? 775 h for z e Si, (6.9) 



-7/7? _ 7a? , for z G & 



Remark 6.4. The asymptotic condition (|6.8I) for X comes from the asymptotic 
condition in the RH problem 15.41 for Y combined with the asymptotic behavior 
(|5.19l) of the Airy functions. The condition (16.81) . however, is not valid uniformly 
as z —> oo for z close to E 2 . More preceisely, (I6.8[) holds uniformly as z — s- oo 
with 



argz e (— 7r + £,— 7r/3 — e)U(— 7r/3 + e,7r/3 — £)U(7r/3 + e,7r — e) 
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for some e > 0, but not for e — 0. This is due to the fact that the functions 1/3, 
2/5 and 2/4 appearing in the second column of the transformation (|6 . 5[) (which 
are dominant solutions of the Airy equation in the respective sectors used in 
(I6.5I0 are not the recessive solutions in the neighboring sectors. 

This non-uniformity is a minor issue that will be resolved more or less auto- 
matically during the transformations of the steepest descent analysis, see also 
[14j[T5]. To fully justify the analysis it suffices to supplement the RH problem 
for X with the following condition that controls the behavior as z — > 00 near 
£2. 

• As z — > 00 we have 



X(z) = 0(l)A(z) 







z ~ n / 2 e 3t °^ T 3~ 



.3/2 



x < 



\0 

z n 



( 











In .3/2 

z -"/2g 3t 0v /tJ J 

\ 



z~ n l 2 e 3t ov / *3 

\ 



.3/2 





2n .3/2 

z -n/2 e 3t 0v /t3~ J 



with O(l) term that is uniform in z. 
Remark 6.5. If we put 



^a 2 3 
e 3t o 











1 














z e S , 

zeSiU S 2 , 

(6.10) 



then it is easy to verify that $ satisfies a RH problem with piecewise constant 
jumps. Hence $ satisfies a differential equation 

±*{z) = C(zMz) 
dz 

with a polynomial coefficient matrix C(z). In addition there are also differen- 
tial equations with respect to to and £3, as well as a difference equation in n. 
We will not elaborate on this here. Instead we are going to make a different 
transformation. 



6.2 Second transformation 

± ^_ z 3/2 

In the next transformation the factors e 3 *ov'3 are removed from the asymp- 
totic condition (|6.8[) in the RH problem for X. The transformation also has the 

TitO 3 

effect of eliminating the (1,3) entry ^pie 31 ? 2 from the jump matrix (|6 . T[) on 
the unbounded contours Cj. 
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Definition 6.6. We define 



/ i ™*3^3\ 

i o i e w 

V(z) = X(z)\o i o 

VO 1 



(6.11) 



for z in the domain bounded by Co U G\ U C*2 (this is the unbounded domain 
containing £2 in its interior, see Figure [5]), 



and 



= X(z) elsewhere, 



\ 





(6.12) 



V(z) = V(z) x < 





f 1 










2n 
g 3t 











(1 










2» 3 










in So, 



(6.13) 



in Si U S2. 



e 3 «oV*a y 


_2 _n_ z 3/2 

e 3 *ov l 3 y 

By straightforward calculations based on the RH problem 16.31 for X and 
the definitions (|6.11l) — (|6.13|) we find that V is the solution of the following RH 
problem. 

RH problem 6.7. The matrix valued function V : C \ T v -> C 3x3 satisfies 

• V is analytic in C \ Ty where Ty = Fx- 

• V-\- = V_ Jy on Ty with 



/ 1 e *o ^ytj 1 1 3 n \ 





; 



v 



J v {z) = < 



4" -3/2 

w 2 e 3 t0v /i? 



(|-±n/3)e t0 3 v^ 



_( 2 Z 3/2_M 3) 
'0 V 3./t^ 







±„ 3/2 



*" -3/2 

w 2 e 3t oV^ 



arg z = 7r/3, 

_ a f 2 3/2 , t_s 3A 

(i - §iV3)e*° l V5* +32) \ 

1 

4 



argz = — 7r/3 or argz = 7r, 



2 3/2 



J 



(6.14) 
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and 



id 



Jv{z) 



(i±iiv^)e 
1 



n I 2 .3/2 . % 3-, \ 

1 (§±±n/3V°^ +32) 0\ 
1 



* G C ± , 



IV 







V 



(6.15) 



'z a 
^(z) = (I+0(l/z))A(z) | z""/ 2 J as : , x. 

z-"/ 2 / 



where A(z) is given by (|6.9p . 

In (|6.14[) we used the convention that arg 
z 3 / 2 = — z|z| 3 / 2 for argz = ir. 

.3/2 



r 3/2 



argz, which means that 



The entries e 3t ov / *3 m the jump matrix (I6.14[) on £2 are oscillatory. 
Later they will be turned into exponentially decaying entries by opening up un- 
bounded lenses around £2- This indeed leads to exponentially decaying entries 
because of the upper triangularity of the jump matrices. 



6.3 Third transformation 

In the third transformation we make use of the minimizer (/Lt*,//^) of the vec- 
tor equilibrium problem. Associated with the measures [i\ and [i\ are the g- 
functions 

g k (z)= f lag(z - s)d»Us), k = l,2, (6.16) 

with appropriately chosen branches of the logarithm. We choose the branches 
in such a way that 171 is defined and analytic in C \ (Si U R~) with g\{x) real 
for x real and x > x* . Then we have the symmetry 



gi(uz) = g 1 (z) + 2?ri/3, 
gi{uj 2 z) = gi(z) - 2ni/3 



, , zeS . (6.17) 



The branches of the logarithm in the definition (|6 . of gi are chosen such that 
g2 is defined and analytic in C \ £2 = So U Si U S2, and g 2 {x) is real for real 
x > with the symmetry 



92(WZ) = 52(2) + 7Tl/3, 

# 2 (w 2 z) = g 2 {z) - iri/3, 



2 it zeS Q . (6.18) 



The conditions (|4. ID — (I4.2[) lead to the following properties of the ^-functions. 
We emphasize that all contours are oriented away from and towards 00, so 
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that for argz = w, we have that gi,-\- (z) is the limiting value from the lower 
half-plane. 

Lemma 6.8. (a) We have, with the constant i from (|4.2[) . 



9i,+( z ) +91 ,-(*) - 52(2) 



3io v^i^ 
2 



z e [0,x*], 



,3/2_^3 +£i 



(6.19) 



■«- = 7r« 3/2 -^-* 8 + ^- 



7TI, Z G [0, W 2 X*] 



(b) PFe ftave 



.92, +(2) +52,-(z) - 51(2) = 0, argz 
92,+ (z) +g 2 -(z) - gi,±(z) = ±ni, argz 



±tt/3, 



(6.20) 



Proof, (a) The equality of the real parts follows from (|4.2[) . since supp(/4) 



Si = U,[o,< 



Both sides of (|6.19|) are real for z £ [0, x*] and so the identity 



(|6.19[) holds on [0,x*]. The identity on the other intervals, then follows from 
the symmetry (|6 . 1 T|) . 

(b) The identities in f|6 . 20[) follow in a similar way from (14.11) . □ 



With the g-functions and the constant I we make the next transformation. 
Definition 6.9. We define 



f e -nt 







1 e -n{gx{z)-l) 







1 


I V(z) 




e n(gi(z)-g 2 (z)) 










V 









g«92(z) 



(6.21) 



It is our next task to state the RH problem for U. It will be convenient 
to write the jump matrices in terms of the two functions tpi and </? 2 defined as 
follows. Recall that £1, £2 an d £3 are defined in (|4.6[) - (|4.8|) . 

Definition 6.10. We define (fx 

1 



\ (Ei UE 2 )^C by 
- &(s))ds, zeSjMO^x*], j = 0,1,2, (6.22) 



2^0 Jujix* 

with integration along a path lying in Sj \ [0,uj j x*] if z G Sj \ [0,u} J x*] for 
j = 0, 1, 2, and tp 2 : C \ {z G C | z 3 G R} ->■ C by 

— 7ri/6, < argz < 7r/3, 



^ 2(Z) = 2^ 



+7ri/6, 7r/3 < argz < 27r/3, 

+7ri/3, 27r/3 < argz < 7r, 

+7ri/6, — 7r/3 < argz < 0, 

— ni/6, — 27r/3 < argz < — 7r/3, 

— 7rz/3, — 7r < argz < — 27r/3, 



(6.23) 



with integration along a path in the sector kir/3 < args < (k + 1)tt/3 if z lies 
in that sector, for k = —3, . . . , 2. 
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The basic properties of ipi and (p2 that connect them with the g-functions 
are collected in the next lemma. 



Lemma 6.11. (a) For z £ Si, we have 

.9i,+0) - 9i-{ z ) = ±2<pi,±(z). 
(b) For z e C \ (Si U S 2 ) we have 



(6.24) 



2g 1 (z)-g 2 (z)-t = 2<p 1 (z)± 

with + in So and — in S\ U S2 ■ 
(c) For z e S 2 , we /lave 



z 3/2 _ ^3 

3i 







in So, 



iri in S±, 
—iri in S2, 



(6.25) 



92,+ ( z ) - 92,-(z) = 1ip 2 ,+{ Z ) ± 



3^0 V^3 



,3/2 



= -2p 2 ,-(*)± 



,3/2 



(6.26) 



3^0 

w/iere we wse + /or argz = 7r/3 and — for argz = — 7r/3 or argz = 7T. 
(d) For z <G S2, we ftawe 



#i,-( 2 ) +32,+ 

f 2^,_ (z) + ^^ 3/2 - |^ 3 , argz = tt/3, 



= < 



y 3/2 t 3 ~3 



3t 



2 . 



arg z = 7r, 



(6.27) 



l 2 ^.-W - ^7H^ 3/2 - - ™> argz = -tt/3. 

(e) FFe have for z £ Si, 

'0, /orze[0,x*], 
2pi i+ (z) + 2<p 2 ,+ (2) = 2ip 2 ,-( z ) + { ™, f° r z € [0,^x*], 

7tz, /or z € [0, cj 2 x*]. 



(6.28) 



and for z G S 2 , 

2<pi i+ (z) + 2^2, + (2) = 2(^i _(z) 



7ri, /or argz = ±7r/3, 
27ri, /or argz = ir, 



(6.29) 
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Proof. By (|4~4l and (ETTo) we have that g[ = F x and g 2 = F 2 , so that by 
(|4.8p we can express the derivatives of the left-hand sides of f|6 . 24[) (|6 . 26[) in 
terms of the ^-functions. This yields 

<?i,+(X> - a'iA z ) = U^+( z ) - (*))> z e s x , 



1 2z X /2 4 



2 fl U*)-fl' 2 (z) = -(&(*) -&(z))± 



3 .2 



Z 



to tovh to 

£0 toV< : 3 



with the same conventions on ± signs as in parts (b) or (c) of the lemma. In 
view of the definitions (|6.22p - (16.23l) we have 



2^-00 = i&(*)- &+i (*))> 3 =1,2. (6.30) 



Using also that £i ; ± = ^2,^ on Si and £3^ = ^2,=f 011 ^2, we then easily check 
that the derivatives of the two sides in (|6.24p - (|6.26l) agree on the respective 
contours or regions. It follows that (|6.24[) - (I6.26[) hold up to a possible constant 
of integration. 

Both sides of (|6.24[) vanish for z — uj^x*, j — 0,1,2, while the equality in 
(|6.25j) for z = uPx*, j — 0,1,2 follows from the identities f|6 . 19[) . Hence the 
equalities in parts (a) and (b) follow. 

To establish the equality in part (c) we examine the behavior as z — > 0. 
Because of the symmetry (|6.18|) we find that there exists a constant 72 such 
that g 2 has the following limits a the origin 

{72, z e s , 

72+™/3, z€Si, (6.31) 
72 - 7ri/3, z e S 2 . 

From (|6.3ip we find the limits of g 2 ,+ {z) — g 2 .-{z) as z — > on the three half- 
rays in £2. These values correspond to the limits of ±2(p 2j ±(z) as z — > that 
we obtain from the definition (|6.23p and part (c) follows. 
Using flSHD]) we find for z £ E 2 that 

/ x, ,s J2ffi0) -9%,-{z) -I, argz = ±7r/3, 

gi,-(z) + g 2 ,+{z) -£= < . (6.32) 

[2.9i -{z) - g 2 -{z) -l-m, argz = n. 

Then (|6 . 27[) follows by letting z in (I6.25P go to the minus-side of £2 , and inserting 
this into (f02|) . 

Because of (|6.30j) we have 

to 
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Using = £2,- on ^i an d £3,+ = £2,- on £ 2 , we then find that the derivatives 
of the two sides of (16.28[) and (|6.29[) agree. To prove the two equalities we again 
examine the behavior as z — > 0. From f|6.31[) and (I6.20j) we find the limits of g\ 
at the origin, 

{27 2 + 7ri/3, < argz < 2tt/3, 
272 + ^ 2 7 r/3<argz< 7 r, 
2 7 2 -vri/3, -2tt/3 < argz < 0, V ; 

272 — 7ri, — 7r < argz < — 27r/3, 

where 72 is the constant from (|6.3ip . Using (|6 . 10[) we also find 

^ = 3 72 . 

Combining (|6.31[) and (16.33)) with (|6.24p and (|6.26l) we obtain the limits of (p% 
and if2 at the origin. 

= ... if fcTT/3 < argz < (fc + IK/3 with k even, 
■ ' I -2tt«/3. if /br/3 < argz < (fc + 1)tt/3 with fc odd, 

< argz < 7r/3 or — 27r/3 < argz < — n/3, 

Inn 2 ,. :>=<"' " < ar § 2 < OT ^ S < < (6 .35) 

^0 I 27r£/3, 2tt/3 < argz < tt, 

[ 2ni/3, —it < argz < — 27r/3. 

The values in (|6.35|) also follow from the definition (|6.23|) of tp 2 - From (|6.34[) 
and (I6.35[) we check that the left- and right-hand sides of (|6.28[) and (|6.29l) have 
the same values as z — > and the equalities (|6.28l) and (I6.29P follow. □ 

Then we have the following. 

RH problem 6.12. U is the solution of the following RH problem. 

• U is analytic in C \ IV, where Tjj = Ty, 
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• U+ = U-Ju on Tjj with 

/ e -2n<pi 1+ (z) 



e -2n V i,-{z) o ( 



J v (z) = < 





V 







(\ 


e 2n Vl (z) q\ 







1 






l) 











2 


a; 2 e -2n V2 ,+ (z) 











fi 


(±±^/3)e 2 







1 










1 

,,,p-2nip2,-(z) 



Z € El, 



z G E 2 , 



± 



(6.36) 

• U(z) = (I + 0(l/z))A(z) as z — » oo, where A(z) is given by (I6.9[) . 

Proof. From (16. 14)) on (J • (0, a; J 'x] and the definition (|6.2ip we obtain by the fact 
that g2 is analytic, that Jjj on [J^- [0, Wjx] takes the form 



e »(si,+ -9i,-) 





3/2_*a._3 



«*)) 



(6.37) 



The (1,2) entry of (|6.37p is equal to 1 on Ei because of (|6.19l) (recall that n 
is even by our assumption at the beginning of subsection 16. 1[ so that the terms 
±iri that are in (I6.19|) on [0,ojx*] and [0,o; 2 a;*] do not matter). The nontrivial 
diagonal entries in (|6.3Z[) reduce on Ei to the ones as given in (|6.36l) on Si 
because of the identities (|6 . 24|) . On \Jj[u^x*,uPx] the diagonal entries in f|6 .37[) 
are equal to 1, since gx t + = g\- is analytic. The (1,2) entry is e 2n(pi because 
of the identity ([BTSSj) . 

For the jump matrix Jjj on £2 we obtain from (|6.14j) . (|6.2ip and the fact 
that e ngi is analytic on S2 



/l 





_ 9 -n(S2,++S2,-± 

ore 



3*0 



"(ff2,+ +32,-± — — 7= 



,3/2-, 



Vo 

(6.38) 

with + for argz = 7r/3 and — for argz = — 7r/3 or argz = n. The nonconstant 
entries in (|6 . 38[) are simplified using (|6.26|) and (|6.27|) which leads to the expres- 
sions given in f|6 . 36[) . Again we use the fact that n is even so that the term ±7ri 
in (|6.27|) on argz = — 7r/3 does not contribute. 
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For z £ (J • we find from (|6.15p and (16.21)) that the jump matrix Jy is 
equal to 

1 

v o 

Here we simplify the (1,2) entry by means of (|6.25[) and thereby obtain the 
jump matrix that is given in (|6.36[) on (J . Cj . 

The asymptotic condition in the RH problem 16.121 follows from the asymp- 
totic condition in the RH problem 16. 71 for V, the transformation (|6.21[) and the 
fact that 

gx{z) =logz + 0(z- 3 ), g 2 (z) = § log 2 + 0{?-^ 2 ) 
as z — > oo. □ 

The jump matrices Jjj in the RH problem 16.121 for U have a good form 
for the opening of lenses around Si and E 2 . Indeed, the expressions <Pj,±(z), 
j = 1,2, that appear in the non-trivial diagonal entries in the jump matrices 
are purely imaginary and the diagonal entries e~ 2nipi - ±lyZ ^ are highly oscillatory 
if n is large. The opening of lenses will turn these entries into exponentially 
decaying off-diagonal entries. 

The other non-constant entries in the jump matrix Jjj have a factor e Vl 
or e 2nVl '- . These are exponentially decaying provided that Retpi < on 
Uj(( w ' 7 x * t^x] U Cj) and Re<^i,_ < on E 2 . The first of these inequalities 
indeed holds, but the second one is not satisfied on the part of E2 in a neigh- 
borhood of the origin. Fortunately, we can remove the (1,3) entry in the jump 
matrix (|6 . 36[) on E 2 by a transformation in a larger domain that is similar to 
the global opening of lenses in [3]. We will do this in the next subsection. The 
transformation will not affect the jump matrices on Si and \^j[u) 3 x*,u) 3 x\. 

6.4 Fourth transformation 

The transformation U i-> T is based on the following property of the set where 
Re<^i < 0. Recall that ipi, see (|6.22[) . is defined and analytic in C \ (Si U E 2 ). 
From {OH) and (EHo) it follows that 

ipi,±(z)=±irinl([z,u j x*]), for z e [0,u) j x*], j = 0,1,2, (6.39) 

so that Keipi.± = on Si. On S 2 we have the identities (|6.29[) with ip2,+ purely 
imaginary. Therefore Re tpi extends to a continuous function on C, which we 
also denote by Retpi. With this understanding we define the domain D as 
follows 

D = {z £ C I Re ipi (z) < 0}. (6.40) 

Since Re</?i is harmonic in C \ (Si U S 2 ), we have that the boundary dD 
consists of a finite number of analytic arcs that start and end on Si U S2 or at 
infinity. 
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5 



Figure 6: Curves Keipi(z) — and the domains where Re<pi(z) < and 
Ke<fi(z) < 0. The figure is for the values to = 1.8, £3 = 0.25. 

From (|6.39[) it follows that ±itpi^± is real and strictly increasing with a 
positive derivative along each interval [OjClPx*] in Si. Then by the Cauchy- 
Riemann equations, it follows that the closed set C \ D contains the intervals 
[0, uj-'x*), j =0,1,2 in its interior. Note that is also an interior point, since 
Hi has a positive density at 0. Near the branch point x* we have 

tp^z) = -c(z - x*f /2 (l + O ((z - x*)" 1 )) asz^i* 

with a positive constant c > 0. This follows from the fact that fi\ vanishes as 
a square root at x*. Thus Re(£>i < immediately to the right of x* . It also 
follows that x* £ dD and dD makes angles ±27r/3 with [0,x*]. The same holds 
at uj j x* because of the rotational symmetry. 
Since 

2ipx{z) = — z 3 + 0{z 3/2 ) as z ^ 00, 

which follows for example from (|6.25|) , the domain D is unbounded and extends 
to infinity in the three directions 

|argz- (2fc- 1)tt/3| < tt/6, k = 0,1,2, 

around the directions ±7r/3, n of £2- The unbounded parts of dD then extend to 
infinity at angles ±7r/6, ±7r/2 and ±57r/6. From what we now know about D it 
follows that 3D has three unbounded parts. One of them intersects the positive 
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real line and tends to infinity at angles ±7r/6. The other two are obtained by 
rotation over angles ±27r/3. 

The analytic arcs in dD that start at the branch points w J x* then necessarily 
end at points on £2. The conclusion is that we have a situation as shown in 
Figure [51 In particular the domain D is connected. 

We also find that Re tp\ takes on a negative minimum on the interval (x* , 00) 
which by (|6.22l) and parts (c) and (d) of Lemma l4~3l is taken at x. Thus x e D, 
and by symmetry also uflx £ D for j — 1,2. We may then also assume that the 
unbounded contours Cj from (|5.13p - (|5.14l) (see also Figure SJ) satisfy 

CjCDHSj, j =0,1, 2, 

and that they extend to infinity at asymptotic directions — tt/3 + 9q and 7r/3 — #0 
(for Co), 7r/3 + 9 and tt — 9 (for C%) and — n + 8 and — ir/3 — 9 (for C2), 
where € (0,tt/6) is some fixed angle. The above choice makes sure that 

Re^i(z)<0, z e Cj, j = 0,1,2, (6.41) 

and, with positive constants c\, C2 > 0, 

Re^i(z) < -ci|z| 3 + c 2 , (6.42) 

for z in the domain bounded by Co U C\ U C2 . 

After these preparations we define the transformation U H> T as follows. 

Definition 6.13. We define T as 

f (z) = U{z) 

for z outside the domain bounded by Co U C\ U C2 , and for z inside this domain, 

,u 



T(z) = [/(z) I I I for , a . rgZ G ( ^' 



and fc = —3, . . . , 2. 

(6.43) 



Then T is defined as 



5 ON 




T(z) = I 1 I T{z) I 1 I . (6.44) 
,0 1, 

The RH problem for T is the following. 
RH problem 6.14. T is the solution of the following RH problem. 
• T is analytic in C \ Tt, where Tt — Tjj, 
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• T + = T-Jt on Tt with 



Mz) = { 




1 




— 2nv2,-(z) 



z e Si 



z G S 2 , 



e U L,E 



una; , w J x 



(6.45) 



as 



• T(z) = (I + 0{l/z))A(z) as z -t oo, where is given by (l6J|) . 

Proo/. From (|Q2]l and (|Q3]) we get that T(z) = U(z) (i + O (e~ Cl ^l 3 )) 
2 — >• oo, which by the asymptotic condition in the RH problem l6.12l for U leads 
to 

T(z) = (I + 0(l/z))A(z) 

as z — > oo. The transformation (|6.44[) does not affect this asymptotic behavior, 
because of the special form (|6.9j) of A(z). 

We next verify the jump matrices (|6.45|1 . We write 



1 1 i— 
a=- + -iV3. 



(6.46) 



For the jump matrix Jf = ( T_ ) T + on Ex we obtain by (|6.36[) and (|6.43l) 



'1 ae 2n ^ 
1 




0> 



'1 ae 2nipi -+ 



o I ./, j 



f e -2n<px,+ l + a + a 0^ 



-2n<pi 






o 



The transformation (|6.44l) has the effect of dividing the (1, 2) entry in Jf by 2 
and we obtain the jump matrix Jt on Si as given in (|6.45l) . Similar calculations 



lead to the jump matrices Jt on (J ■[* 



U) J x , uj j x 



and on [j • C. 



± 

j "j 
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For the jump matrix on S2 we obtain by (|6.36l) and (I6.43[) 




(6.47) 



By (RT2g|> we have e 2«(vi,--V2,+) _ e 2n<pi,+ on j] 2 , and so the (1,2) entry in 
(|6.47p reduces to 

{-a - aco 2 )e 2nVl '+ = 0, 

because of the value (|6.46p we have for a. Thus the (1,2) entry in (|6.47l) 
vanishes. The transformation (|6.44l) then gives Jt — Jf and it is equal to the 
jump matrix (|6.45|) on S2- □ 

Since Ret^i < on [_]-{{uj^x*,ijj^x\ U C^) we see from (|6.45|) that the jump 
matrices Jt on these parts of Tt tend to the identity matrix as n — > 00. What 
remains are oscillatory jump matrices on Si and S2. 



6.5 Fifth transformation 

We open up lenses around the three intervals in Si and around £2. We use 

L 1= L+ULi (6.48) 

to denote the lens around Si where L\ (L±) is the part that lies on the + side 
(— side) of Si. 

The boundary dL\ of L\ meets Si at the endpoints uj^x* , j — 0, 1, 2 but not 
at 0. The boundary meets S 2 at points at a positive distance £1 > from the 
origin, see Figure [71 and we choose L\ such that dL\ \ {ui^x*} is contained in 
the region where Re (fx > 0, see also Figures [6] and [Jj 

Similarly we write 

L 2 = L+ U Lz (6.49) 

for the lens around £2. The boundary dL 2 lies in the region where Re if 2 > 
and they are bounded by infinite rays that meet Si at points at a positive 
distance 62 > from the origin. The rays have an asymptotic angles ±7r/3 ± e, 
7r ± e for some small e > 0, see also Figure [7J We can indeed open the lense 
around S2 this way, since n* 2 has a positive density, also at 0. 

Note that the lenses L\ and L2 have non-empty intersection. More precisely, 
L+ n ^ and D L+ ^ 0_ 

Definition 6.15. We define 

S(z) = T(z) 
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Figure 7: Lenses L\ and L 2 around Si and £ 2 and the contours E5 in the RH 
problem for S. 



for z outside of the lenses L\ and L 2 - For z inside the lenses we define 



S{z) = T(z) x < 



1 

Te -2n Vl (z) 1 J 

i / 
/l o\ 

1 , 

1 o\ 

- Fe -2n Vl ( z ) x Q 

y w T e -2n( ¥ >i(z)+ip 2 (z)) _|_ u; ± e -2n ¥ ) 2 (z) j 



z e Lf \L 



2- 



(6.50) 



where w+ = w and w = w 1 = u 2 . 



Then S 1 satisfies a RH problem on a union of contours T$ that consists of Tt 
and the lips dL\ and <9L2 of the lenses. We continue to follow the convention 
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that the contours are oriented away from the origin, and towards infinity. For 
dh\ and dh-i this means that each part in dL\ is oriented towards a branch 
point uj 3 x* , and each part in dL*2 is oriented towards infinity, as indicated in 
Figure [JJ 

RH problem 6.16. • S is analytic in C \ Ts, 

• 5+ = S— Js on Ts with J$ given as follows on IV, 

r / o i o> 

-1 | , z G Si, 




zes 2 , 



(6.51) 



U! J X ,U} J X\ 



and on the lips of the lenses by 

1 

e-^iW 1 I , 





1 o x 
1 

Q cjTe -2« V2 ( z ) ^ 
1 

± e -2n(tpi(z)+tp 2 (z)) q ^ 

1 N 

1 



z e 3Li \ £ 2 , 
2 G dL±\L u 

z e dLf n L 2 , 
z e dL* n Li. 

(6.52) 



• 5(z) = (1+ 0(l/z))A(z) as z -> oo where A(z) is given by (|6T9|) . 

Proof. The jump matrices J5 arise from combining the jump condition (|6.45l) 
for T with the definition (|6.50p . As an example we show how the jump matrix 
(|6.52p on dLf n L 2 arises. 
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To compute the jump matrix on dL^ n L 2 we need the definition (|6.50[) of 
S in the parts L\ n L 2 (which is on the —side of dh\ n L 2 ) and \ Li (which 
is on the right). Since T has no jump we find from (|6.50l) that, 

1 °\ 1 A ° N 

J, A z) = I e -2n<PiM 1 I [ I 





since 1 + cu + ui 2 = 0. This then easily reduces to the jump matrix Jg as given 
in (|6.52p on dL-^ n L 2 . The other jumps follow in a similar way. 

The asymptotic condition for S outside of the lens L 2 follows trivially from 
the asymptotic condition in the RH problem [6T4] for T. Inside L2, we note that 
by (|6.26[) we have 

2 

tp 2 (z) = ± — z 3/2 + O(logz) as z -> 00, z e L 2 , 



where the ± signs are such that Reip 2 (z) > for z large enough in L 2 . Then it 
follows from (|6.50|) that the asymptotic condition is also valid in the lens L 2 . □ 

We are now in a situation where all jump matrices Js tend to the identity 
matrix as n — > 00, except for those on £1 and £2, which are constant. We 
already saw this for the jumps on \J- (oj^x* , uijx]UCj ) , see (|6.5ip , since Re ipi < 
on these parts of the contour. Also L\ is contained in the region where Re ipi > 
0, so that the entries e ~ 2 " Vl ( 2 ) that appear in (|6.52l) are indeed decaying as 
n — » 00. 

Regarding ip 2 , we note that by (|6 . and (|6.26p we have that <p 2 ,±(z) is 
purely imaginary for z S £2, with 



W2,+{z) = -iif2,-(z) 



^([0,z]) + — ^=|z| 3 / 2 



'-tt/6, argz = ±7r/3, (6.53) 



3to^/t3 [k/S, argz = 7r. 

Thus i(p 2j + and —iLp 2> - are strictly increasing along each of the rays in £2. Then 
by the Cauchy-Ricmann equations, we have that Re 922 > both to the left and 
to the right of £2. We may (and do) assume that the lips of the lens L 2 is 
contained in the region where Re (p 2 > and it indeed follows that the jump 
matrices (|6.52p tend to the identity matrix as n — > 00. 

Remark 6.17. The asymptotic condition in the RH problem 16.161 for S is 
valid uniformly as z — > 00 in any direction in the complex plane. This is in 
contrast to what happens in the RH problems for X, V, U and T, where the 
asymptotic condition is not valid uniformly in the directions argz = ±7r/3, 7t 
that correspond to £2, see Remark 16.41 
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The term e ~ 2nV2 ^ that appears in the transformation (I6.50[) for z 6 L^\L\ 
is small for z away from £2, but becomes of order 1 as z approaches £2- On 
£2 a certain cancellation takes place, which results in the asymptotic condition 
being valid uniformly. One can establish this rigorously by looking back at the 
transformations Y i-> X V i-> U M>Tn-S' near £2, but we will not go into 
that here. 

If we do not want to rely on this detailed analysis, then we can only guarantee 
at this stage that the asymptotic condition in the RH problem l6.16l for S is valid 
uniformly as z — > 00 with 

argz 6 (— 7r + e, — n/3 — e) U (— 7r/3 + e,7r/3 — e) U (n/3 + e,it — e) 

and that 

S(z) = 0(l)A(z) uniformly as z -> 00. (6.54) 

The latter condition is a consequence of (|6. 10|) in Remark 16.41 and the transfor- 
mations leading from X to S. 



6.6 Global parametrix 

The global parametrix M is a solution of the following model RH problem, which 
we obtain from the RH problcm l6.16l for S by dropping the jump matrices that 
tend to the identity matrix as n — > 00. 

RH problem 6.18. • M is analytic in C \ Tm where Ym = £1 U £2, 
• M + = M-Jm on V m with 



J M {z) 




z e £1, 



z e £ 2 , 



• M(z) = (1+ 0{\/z))A(z) as z ->• 00, where A(z) is given by (|6.9p . 

• M(z) = O ((z - ujix*)- 1 / 4 ) as z -> for j = 0, 1, 2. 

• M(z) remains bounded as z — > 0. 

Note that the product of the jump matrices around is the identity matrix 
and so we can indeed require that M remains bounded near the origin. Also 
observe that the asymptotic condition is compatible with the jumps on the 
unbounded contours, since 

A+=A-J M on £ 2 , (6.55) 
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which may be checked directly from the definition (I6.9P of A(z). 

To construct M we use a meromorphic differential f2 on the Riemann surface 
1Z, which we specify by its poles and residues, see also Q]3 [351 [25]. We require 
that f2 has simple poles at the branch points ui^x* , j — 0, 1, 2 and at 002 (the 
point at infinity that is common to the second and third sheets) with 

Res(f2,c^'a;*) = -~, Res(Q, oo 2 ) = § 

and ft is holomorphic elsewhere. The residues add up to 0, and therefore the 
meromorphic differential exists and it is also unique since the genus is zero. 
We use 00 1 as base point for integration of and define 



uj(z)= / O, zeTZj, j = 1,2,3, (6.56) 

«/ oci 

where the path of integration is chosen according to the following rules: 

• The path for u\(z) stays on the first sheet. 

• The path for 112(2) starts on the first sheet and passes once through £1 to 
go to the second sheet and then stays on the second sheet. 

• The path for 113(2) starts on the first sheet, passes once through £1 to 
go to the second sheet, and then passes once through £ 2 to go the third 
sheet and then stays on the third sheet. 

• All passages from one sheet to the next go via the — side on the upper 
sheet to the +-side on the lower sheet. 

Then the functions uj are well-defined, u\ is analytic on C\£i, u 2 is analytic 
on C \ (£1 U £2), and u 3 is analytic on C \ £2. The functions satisfy 

m 3j+ = w 3i _ ui l _=ii2,+, "1,+ = «2,- ± 7ri, on £1, 
ux t+ = ui-, u 2 -=u 3 . + , u 2 ,+ = u 3 - ± ni, on £ 2 . 



We put 

and then we have 

(vi v 2 v 3 ) = (vi v 2 v 3 )_Jm on£iU£ 2 , 



(6.57) 



Since v x (z) = 1 + 0(1/ z), v 2 (z) = 0(z- 3 / 4 ), v 3 (z) = 0( Z - 3 / i ) as z -> 00, the 
vector (v\, v 2 ,v 3 ) satisfies the conditions for the first row of M. 

The vector space of holomorphic functions on 1Z \ {002} with at most a 
double pole at 00 2 has dimension 3. Let /« = l,/( 2 ),/( 3 ) be a basis of this 
vector space. We use to denote the restriction of to the sheet IZj. Then 
it is easy to see that 

(vi v 2 v 3 

,, f(2) „ f (2) „ f (2) 
VlJl v 2j 2 W3/3 
A3) „ f (3) ,, f (3) 
VlJi V 2 J 2 U3/3 
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is denned and analytic in C\ Tm with jump B + = B_Jm- Also B{z) = 0{z 1 ^ i ) 
as z — > oo. It then easily follows that dct B = const and the constant is non-zero, 
since the functions /W = 1, f( 2 \ Z*- 3 ' are independent. 

We already noted that A + = A-Jm on S 2 - Since A and B have the same 
jumps on £2, it follows that BA~ X is analytic in C \ Si, and therefore has a 
Laurent expansion at infinity. Since both B(z) = Oiz 1 / 4 ) and A(z) = 0{z x ^) 
asz-> 00, we have 

B(z)A- 1 (z) = C + Oiz- 1 ) asz^oo 

for some constant matrix C. Since dct£> = const 7^ and det A = 1, we see 
that detC 7^ 0, and so C is invertible. Then 

MO) = C^Biz) 

satisfies 

M(z) = (I + 0(z- 1 ))A(z) asz^oo. 

It also satisfies the jump condition M + = M-Jm, it is bounded at and it has 
at most fourth root singularities at the branch points. Thus M is the global 
parametrix we are looking for. 

Remark 6.19. From the above construction it follows that 

M M (z) = vi(z) = e" l(z) , «GC\Ei, 
which in particular means that M\ } i(z) ^ for z € C \ Si. 

6.7 Local parametrices 

The local parametrix P is defined in the disks 

D{oj j x*,S) = {zeC \ \z-uj j x* \ < 6} 
where 6 > is taken sufficiently small. 

RH problem 6.20. • P is continuous on ((J . D^x* , 5))\Ts and analytic 
in QJ.D(wtar,8))\r s , 

• P+= P-J S on T s n Uj D(u>ix*,5), where J s is as in ([63T ]l -([632 ]) 

• P(z) matches with the global parametrix M(z) in the sense that 

P(z) = M(z) (I + Oin- 1 )) asn^oo, (6.58) 
uniformly for z G (J^. dD(ujjX* , (5). 

In the noncritical case to < to,crit the density of /ii vanishes as a square root 
at lo j x* , for j = 0, 1,2, which means that we can build the local parametrix 
P out of Airy functions in a small disk around each of these endpoints. The 
construction is done in a standard way and we do not give details here. 
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Figure 8: Contour Fr for the RH problem for R. 



6.8 Sixth transformation 

Now we make the final transformation. We choose a small disk D{u)i ,x* ,6) 
around each of the branch points, and define 

Definition 6.21. We define 



Then R is defined and analytic outside of the union of Ts with the circles 
around the branch points. Since the jump matrices of S and M agree on S1US2, 
we see that R has an analytic continuation across £2 and across the part of Si 
that is outside of the three disks. Similarly, the jump matrices on S and P agree 
inside the disks, and so R also has analytic continuous inside the three disks. It 
follows that R is analytic in C \ Tr, where Tr is the contour shown in Figure 
[H It consists of the part of Eg \ (£1 U £2) that is outside of the three disks, 
together with the three circles around the branch points. 

The RH problem for R is thus as follows. 




S(z)P(z)- 
S(z)M(z) 



-1 



1 



in the disks around the endpoints uj j x* , j — 0, 1, 2, 
outside of the disks. 



(6.59) 
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RH problem 6.22. • R is analytic in C \ I\r, 

• R + = J?_ Jr on where 

J Af(z) _1 P(z) for z on the circles, 

H ^ ~ \m(z)- 1 J s (z)M(z) elsewhere on T R . 

• R(z) = I + 0(z _1 ) as z -> oo. 

Remark 6.23. Following Remarks 16.41 and 16.171 we note that the asymptotic 
condition in the RH problem 16.221 is valid uniformly as z — ¥ oo with 

argz 6 (— 7r + £,— 7r/3 — e)U(— 7r/3 + e,7r/3 — £)U(7r/3 + e,7r — e) 

for any e > 0. From Remark 16.171 and the transformation (|6.59l) it also follows 
that R{z) = O(l) as z — y oo uniformly. Since the jump matrix Jr{z) for z G Y r 
tends to the identity matrix as z — > oo, see also (|6.62p below, it is then easy 
to show that the asymptotic condition R(z) — I + 0(z^ r ) as z — > oo is in fact 
valid uniformly. 

All jump matrices Jr in the RH problem 16.221 tend to the identity matrix 
as n — > oo. Indeed, because of the matching condition (|6.58l) we have 

J R {z) =I + 0{n- 1 ) 1 for z € dD(u j x*,6) (6.61) 

while the jump matrices on the remaining parts of Tr are exponentially close 
to 

J R (z) =1 + 0(e- c "l z l 3 ), elsewhere on T R . (6.62) 

We conclude two things from (I6.61[) - (|6.62p . The first thing is that the RH 
problem !6.22l has a solution if n is large enough. This indeed follows from (|6.61l) 
(|6.62l) since for jump matrices Jr close enough to the identity matrix the RH 
problem 16.221 has a unique solution which can be written down as a Neumann 
series. Since the transformations 

Y^X^V^U^T^S^R (6.63) 

are invertiblc it then also follows that the solution Y to the RH problem 15.41 
uniquely exists. So in particular the (1, 1) entry exists, which is the orthogonal 
polynomial P n , n . 

The second thing we obtain from (|6.61[) — (|6.62[) is that the solution R of 
the RH problem 16.221 not only exists but is also close to the identity matrix as 
n — > oo. There is an estimate 



R(z) =7 + as n -> oo, (6.64) 
\n(l + \z\)J 

uniformly for z € C \ Tr. The estimate (|6.64l) completes the steepest descent 
analysis of the RH problem. 
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6.9 Proof of Lemma 16.11 

The proof of Lemma 16.11 now follows by unravelling the transformations (|6.G3[) 
and then use (|6.64[) to see the effect on the orthogonal polynomial, since 

P n , n (z)=Y ltl {z), 

see (j5^5| . 

Proof. From the definitions ([BTTllHpTSl) . (|OT|) . (l6Ti5|) - (|Q4l we 

easily see that 



P n , n (z) = X hl (z) = V ltl (z) = U ltl (z)e n9 ^ 
= T M (z)e" 9l(z) , z£C\E v 



(6.65) 



Using ff6~65|) and (j6~50)) we find 

P n , n {z) = S hl {z)e na ^ z \ zeC\Lt. (6.66) 

For z outside of the disks 'x* ,6) we have S = RM by (|6.59p and so by 
(EH 

Sx,i(z) = (1 + 0(l/n))M hl (z) + 0(l/n) 

which, since Mi ! i(z) is analytic with no zeros in C \ Ei, see Remark l6.19[ can 
also be written as 

$1,1(2:) = (1 + 0(l/n))M ltl (z), z e C \ (Li U (J D(wV, 5)), 

and the O-term is uniform. Inserting this into (|6.66[) we obtain (|6.ip uniformly 
for z EC\ {L 1 L)\J j D(uj j x* 7 6)). The lense Li around Ei and the disks D^X'S) 
around the branch points can be made as small as we like. It follows that (|6.ip 
holds uniformly for z in compact subsets of C\ Si and the lemma is proved. □ 

6.10 Proofs of Theorems EH and 1231 

Proof. For to, t% > with to < to^crit we proved that the orthogonal polynomial 
P n ,n exists for large enough n and satisfies (|6.ip . Since Mi : i(z) ^ for z E C \ 
Ei, the zeros of P n . n accumulate on Ei as n — > oo, and this proves Theorem l2.3l 
From (|6.ip and (|6.16l) it follows that, uniformly for z € C \ Ei, 

lim -log|P„,„(z)| = Regi(z) = f log \z - s|d/^(s). 

n— ¥oo n J 

Standard arguments from logarithmic potential theory, see e.g. Theorem 
III. 4.1] then imply that /i| is the weak limit of the normalized zero count- 
ing measures of the polynomials P n ,n as n — > oo. By construction, /i| is the 
first component of the minimizer of the vector equilibrium problem, see also 
Lemma [131 and Theorem I2J5I follows. □ 
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